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DISTRIBUTIONS IN RELATIVISTIC HEAVY- ION COLLISIONS 



By 
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Principal Investigator: Govind S. Khandelwal^ 

ABSTRACT 

An optical model description of energy and momentum transfer in 
relativistic heavy-ion collisions, based upon composite particle multi- 
ple scattering theory, is presented. Transverse and longitudinal 
momentum transfers to the projectile are shown to arise from the real 
and absorptive part of the optical potential, respectively. Compari- 
sons of fragment momentum distribution observables with experiments are 
made and trends outlined based on our knowledge of the underlying 
nucleon-nucleon interaction. Corrections to the above calculations are 
discussed. Finally, use of the model as a tool for estimating colli- 
sion impact parameters is indicated. 
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CHAPTER I 


INTRODUCTION 

With the advent of relativistic heavy ion beams at Berkeley, 

Saclay and Dubna, the experimental situation in high-energy nuclear 
interactions has improved dramatically over the previous era of cosmic - 
ray heavy-ion physics. Sophisticated experiments have challenged 
theorists to come up with new theoretical tools and insights to under- 
stand the new features apparent in relativistic heavy-ion collisions. 
Traditional nuclear physics has been primarily about the nature of 
nuclear matter at or near equilibrium. With heavy- ion beams, the 
possibility of compressing nuclear matter to two or three times the 
normal density and heating to temperatures - 100 MeV has opened up. It 
is appropriate at this point then to review the salient features of 
high energy heavy- ion collisions beginning with cosmic -ray heavy ion 
physics . 

Ever since the discovery of Z > 2 (where Z is the nuclear charge) 
components in the primary cosmic radiation by Freier et al.l>2, which 
fulfilled a prediction by Alfven-*, the subject of high energy interac- 
tions between nuclei has been of fundamental interest. The goal of 
these pioneering studies focussed primarily on interaction mean free 
paths and reaction cross-sections. The production of nuclear fragments 
and determination of their isotopic composition was intensely studied 
in order to infer from these data the conditions of their origin, 
possible acceleration mechanisms and subsequent propagation. These 
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aspects of cosmic-ray heavy ion physics are thoroughly reviewed by 
Shapiro and Silberberg,^ Waddington^ and in the classic monograph by 
Powell. ^ 

The qualitative classifications of the nuclear interactions in 
cosmic rays were performed first by Bradt and Peters. ^ The concepts of 
" Peripheral" (large impact parameter) and "Central" (small impact 
parameter) collisions were introduced by these authors. In a peripher- 
al collision, part of the nucleus overlapping the target is sheared off 
while the remaining fragment proceeds at near the beam velocity. Both 
projectile and target fragmentation may be described as peripheral 
processes. In a central collision, both the nuclei are destroyed, 
involving high levels of excitation and the emission of large numbers 
of secondary fragments. Nucleons, light fragments and pions are copi- 
ously produced in central collisions. 

Experimental studies of cosmic rays revealed many important fea- 
tures, in spite of low intensities and uncertainties in charge, mass 
and energy determinations. Bradt and Peter s^ analyzed the reaction 
cross sections using a semi-empirical "black sphere" expression 

Vac ‘ ' <4 /3 S *T /3 -S> 2 (I D 

with Ap, Ax the mass numbers of beam and target, and S an overlap 
parameter representing the diffuseness and partial transparency of 
nuclear surfaces. With fixed values of r Q and 6, reaction cross-sec- 
tions could be reasonably predicted with equation (1.1). These cross- 
sections were assumed to be energy- independent for bombarding energies 
from .1 A GeV to 30 A GeV. The average number of fragments produced by 
the fragmentation of a certain projectile was also found to be nearly 
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energy- independent (to within -20%). The Bradt-Peters "black sphere" 
model was later refined and a "grey sphere" model was proposed® to 
account for the reduction of geometric cross-section due to the trans- 
parency of nuclear surfaces. Optical model calculations of reaction 
cross-sections were also undertaken successfully by various authors^ » 
and satisfactory agreement was found with experiment. 

For cosmic rays, the reaction products studied were mainly mesons 
and nucleonsll'15 # The experimental goal was to understand nucleus- 
nucleus interaction as a superposition of independent nucleon-nucleon, 
nucleon-nucleus or alpha-alpha collisions. Although successful in 
achieving a broad understanding of such collisions, precise knowledge 
of high energy interactions of nuclei could not be gained from these 
studies due to low statistics and lack of control over experimental 
conditions . 

The first laboratory acceleration of relativistic heavy- ions was 
accomplished in 1971 at the Princeton Particle Accelerator (PPA) , 
shortly followed by the Berkeley Bevatron^®. Acceleration of alpha 
particles began in Dubna 1970. Similar proposals were made at CERN for 
heavy- ion experiments^?. With the closing of PPA in 1972, Berkeley 
BEVALAC became the only high energy heavy ion accelerator in the U.S.A. 
Proposal for a Relativistic Heavy Ion Collider (RHIC) to be built at 
Brookhaven National Laboratory in the 1990' s is underway. Projects 
that are also in progress are GANIL in France, Numatron in Japan, 
Nuklotron in the U.S.S.R. and GSI in Darmstadt, West Germany. 

The BEVALAC, proposed by Ghiorso^®, employs an 8.5 A MeV heavy ion 
linear accelerator, the SUPERHILAC, to inject the ions into the 
Bevatron which continues the acceleration of these ions to a maximum 
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energy of 2.6 A GeV. High intensity Uranium beams have recently been 
accelerated at the BEVALAC. 

Experimental techniques at high energy heavy ion accelerators 
combine tools from both traditional low energy nuclear physics as well 
as high energy particle physics because of the wide range of charge, 
mass and energy. These techniques range from particle identification 
by AE-E and time of flight over 4* steradians for target fragmentation 
to high energy techniques such as magnetic spectrometers, measurements 
of dE/dx, rigidity, Time of Flight (TOF) for slow projectile parti- 
cles. For relativistic particles Cerenkov radiation as well as the use 
of streamer chamber or other multiple track detectors are used. Detec- 
tion by emulsions, plastics and AgCl monocrystals are also used because 
of their wide range of sensitivities, versatility and small demand for 
beam time. The availability of wide ranges of charge in heavy ion 
experiments pose difficulties for charge identification by dE/dx*Z^ 
f (/9) so that additional capabilities must be incorporated into the 
system. Excellent reviews of these aspects of detector development are 
the Heavy Ion Study Proceedings published by Lawrence Berkeley Labora- 
tory (and GSI) every two years, where recent information on both exper- 
imental and theoretical aspects can be found. 

In this work, we shall examine one aspect of relativistic heavy- 
ion collisions in detail. We shall formulate a theoretical framework 
to describe how momentum and energy are transferred to relativistic 
heavy ions. An optical model description of momentum and energy trans- 
fer based on the multiple scattering theory of nucleus -nucleus colli- 
sions will be presented. 
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Aside from the fundamental importance of the topic itself, the 
necessity for understanding these processes (energy* momentum transfer) 
arose due to the increasing sophistication of experiments. Single- 
particle inclusive experiments of the type Projectile + Target -► Pro- 
jectile Fragment [+ other products] form the basis of our knowledge on 
heavy ion projectile fragmentation. It was observed experimentally^ 
in the fragmentation of ^ C and ^0 (2.1 A GeV) beams on targets 
ranging from H through Pb that the isotope production cross-sections 
were factorizable into 

F 

a FBT “ 7 B 7 T (1.2) 

where ^fbT the cross-section for producing the fragment; 7 ^ and 7 ^ 
are two terms that depend on the beam- fragment and target respectively. 
This suggests that the momentum of fragments in the projectile rest 
frame should also exhibit independence of target structure and beam 
energy. Indeed, in the rest frame of the projectile, the longitudinal 
momentum distributions of fragments show a statistical Gaussian depen- 
dence. Irrespective of projectile, beam energy (> 1.05 A GeV) and 
target nucleus, the longitudinal momentum P| distributions for all 
fragments from ^ C and projectiles are characterized by Gaussian 
shape with rms width a(P||) * 50 to 200 MeV/c. These longitudinal 
spectra are also downshifted by <P |j> * -20 to -130 MeV/c from the 
beam, showing that the mean velocities of these fragments are less than 
that of the beam. The rms widths o(P||) an< * a (Pj.) °f the longitudinal 
and transverse momentum spectra are found to be equal to within - 10 % , 
consistent with isotropic production of fragments in a frame moving at 
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a velocity less than the beam velocity. The widths a ( P | ) and a(P^) , 
independent of target mass and beam energy, depend on the masses of the 
beam and the fragment. A parabolic shape reproduces the general trend 
of the data 

a(P||) - 2o q [x( 1 - x)] 17 ; (1.3) 

A P 

where Ap t Ap are the beam and fragment masses respectively, and o Q is 
experimentally extracted from the data or predicted by theory. 

The parabolic shape of the widths <7(P|j) has heen the subject of 
considerable attention of theorists^- 22 w ho explained this dependence 
using conservation of momentum. According to these theories, the 
Gaussian momentum distributions can be understood by treating the 
fragmenting nucleus as a Fermi gas and assuming (i) Momentum conserva- 
tion, (ii) No correlation among nucleons in the parent nucleus, and 
(iii) Neglecting anti-symmetrization of the single particle states. 
Within these assumptions, one would predict (1.3) where a 0 is related 
to the Fermi momentum of the projectile (Pp) via 



and also that (1.4) 

a(P|j) - a(P ± ) 

explaining the isotropic production of the fragments. These insights 
form the basis of much of the analyses of experimental data. 

The simple Fermi gas picture of projectile fragmentation has been 
questioned by many authors. Nuclear structure and binding-energy 
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effects in determining the momentum widths has been pointed out-^. 
Neglect of anti-symmetrization of single-particle states has been 
questioned by Bertsch^ who showed that the momentum widths would be 
reduced from the free Fermi gas value due to correlations. Murphy^ 
has investigated the constraints of phase space on treating the projec- 
tile as a Fermi gas which emits fragments that are also Fermi gases. 
Recent experimental data^6 on the transverse momentum widths of Lantha- 
num fragments also indicate the inadequacy of the simple Fermi gas 
picture. Indeed, an alternate formulation^? of projectile fragmenta- 
tion abandons the Fermi gas picture and attempts to explain widths in 
terms of neutron separation energies. It is clear that one needs to 
bring in sophisticated theoretical insights to address structure and 
binding energy effects, correlations, phase space constraints etc. 

In this work, we attempt to address important but as yet unre- 
solved questions, i.e. how are momentum and energy transferred to heavy 
ions at high energies and how does energy-momentum transfer affect the 
fragment momentum and energy spectra, fragmentation cross-sections and 
their angular distributions. 

To address the above questions, an optical model description of 
momentum and energy transfer between relativistic heavy ion collisions 
within the multiple scattering theory framework will be presented. The 
inputs into our calculations will be the well-known nuclear ground 
state densities and NN (nucleon-nucleon) two-body transition ampli- 
tudes. The energy -dependence of the two-body amplitude will enable us 
to address the above questions over the entire energy range from low to 
intermediate and extremely high bombarding energies for any projectile - 
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target combination . New insights gained from this work, it is hoped, 
will unify a host of data on heavy ion (peripheral) fragmentation. 

Previously, excitation energy calculations have been undertaken by 
Hufner^S e t al. within the Glauber theory framework. Fricke^ has 
calculated the excitation energy of "anomalons" detected in heavy ion 
experiments using an impulsive excitation picture. Phenomenological-^ 
and semi -empirical estimates of the excitation energy in heavy ion 
collisions also exist (these are not based on two -body interaction 
parameters) . 

The new feature of this work is the introduction of a complex 
momentum transfer vector which results from the use of a complex two- 
body transition amplitude that satisfies unitarity and is used to 
analyze experiments. The real (imaginary) part of the two-body inter- 
action, folded with the appropriate densities for nucleon-nucleus or 
nucleus -nucleus scattering gives rise to the real (imaginary) part of 
the optical potential. The imaginary part accounts for inelastic 
scatterings as well as true absorption, and as will be shown, gives 
rise to longitudinal momentum transfer. The reaction cross-section in 
Glauber theory, for example, is obtained from the imaginary part of the 
phase shift computed from the imaginary part of the optical potential. 
The real part of the complex momentum transfer vector represents trans- 
verse momentum transfer due to elastic scattering at high energy. In a 
fragmentation experiment, the projectile (or target) fragments are 
detected. It has been found that these fragments emerge with less than 
the beam velocity i.e. the longitudinal momentum spectra of these 
fragments are "downshifted" by -20 Mev/c to -130 Mev/c depending on the 
fragment. This "momentum downshift" is naturally explained in this 
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work as due to longitudinal momentum transfer to the projectile by the 
target. This arises due to inelastic scatterings that occur as the 
projectile traverses the target. The transverse momentum transfer due 
to elastic scatterings, is separately obtained from the real part of 
the interaction. Calculation of momentum transfer allows us to calcu- 
late the means and widths of transverse (and longitudinal) momentum 
spectra for fragments for any projectile-target combination, using as 
inputs nuclear ground state densities, NN transition amplitude uti- 
lizing the currently available theories of projectile fragmentation. 

Since the two-body interaction that is used is energy dependent, 
the validity of these insights can be tested at extremely relativistic 
energies (> 10 A Gev) as well as at lower (< 1 A GeV) bombarding 
energies, assuming the validity of the underlying model at these two 
extremes. Silicon beams have now been accelerated to energies “14.5 
A GeV (at Brookhaven) , Oxygen and Sulphur beams to 60 A Gev and 200 A 
GeV at CERN. The momentum spectra of fragments look remarkably similar 
to those at 2.1 A GeV. Lower (< 1 A GeV) energy data on these same 
fragments is harder to come by; the paucity of available data prevents 
a systematic study of the reaction mechanism and its evolution as a 
function of bombarding energy. With the methods described herein, one 
can now perform a theoretical calculation within this model. Such 
calculations have been included in this work for possible future com- 
parisons . 

Another often studied question in nucleon-nucleus as well as 
nucleus -nucleus collisions is the validity of the constant velocity 
assumption that is frequently made. Within this model, we have a 
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computationally tractable scheme for addressing this question. Longi- 
tudinal momentum transfer to the beam nucleus inherently challenges the 
notion of constant velocity. The magnitude of the momentum transfer is 
responsible for change in velocity. We have developed a scheme that 
takes into account corrections to the constant velocity assumption used 
in our formalism. Another correction we address is the possible Cou- 
lomb effects as the electromagnetic fields generated at these energies 
are substantial enough to require such an analysis. 

We have also derived a theory of energy transfer (based on multi- 
ple scattering theory) along the lines of our theory of momentum trans- 
fer. Experiments are seldom able to measure the excitation energy of 
fragments; it is rather inferred or extracted from observed cross- 
sections (using models such as the "Ablation-Abrasion” model). Our 
calculation of excitation energy for projectiles and projectile frag- 
ments introduces a comprehensive framework for performing these energy 
transfer calculations using NN interaction parameters as inputs. No 
comparisons could be made with experimental data because of the absence 
of the latter. However, sophisticated experiments in the future may 
change this situation. 

The remainder of this work is organized as follows: In Chapter 

II, the multiple scattering theory of nucleus -nucleus collisions is 
reviewed and the optical model discussed within this context. In 
Chapter III, the formulation of momentum and energy transfer in 
nucleus -nucleus collisions is made within the multiple scattering 
theory framework. Chapter IV contains numerical results on momentum 
transfer (longitudinal and transverse), momentum downshifts and momen- 
tum widths, preceded by discussions on the theoretical understanding 
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and experimental facts on these topics. Corrections are discussed in 
Chapter V. The first is the "deceleration correction" arising from 
the transfer modification of the constant velocity assumption at high, 
intermediate and low bombarding energies . The second is the more 
familiar Coulomb correction. We conclude by discussing the major new 
insights gained in this work as well as indicating possible directions 


for future research. 



CHAPTER II 


MULTIPLE SCATTERING THEORY OF NUCLEUS -NUCLEUS COLLISIONS 
Multiple scattering theory provides a reliable and fundamentally 
correct description of hadron-nucleus as well as nucleus -nucleus 
collisions at intermediate (- 100 A MeV to 500 A MeV) and high (>500 
A MeV) incident energies . 31-39 j n these theories, the complex many- 
body problem of Ap projectile nucleons interacting with A<j target 
constituents is formulated in terms of two-body interactions. The 
success of these theories is well documented and is reviewed quite 
frequently in the literature. 

Our goal in this chapter is to review multiple scattering theory 
in order to formulate the problem of momentum transfer and excitation- 
energy deposition in nucleus -nucleus collisions. We shall review the 
main results of the multiple scattering theory of nucleus -nucleus 
collisions. Following Wilson, 39,40 a set Q f coupled equations relating 
all entrance channels to all exit channels will be derived. An optical 
potential V 0 p t will be extracted under certain approximations. 

II. a Review of Multiple Scattering Theory 
We shall review multiple scattering theory of nucleus -nucleus 
collisions to find an expression for the multiple scattering series. 
Collision of a composite projectile (mass number Ap) with a composite 
target (mass number Aq0 will be considered. The formulation of this 
problem and its solution can be found in references 39-42. The 
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Hamiltonian for the combined system of N - (Ap + A-p) nucleons can be 
written as 




H - 2 T. + 2 


. V. .+ 2 T + 2 V a 

j-1 J i<j lj a-1 ° <*<P ^ 


+ 2 V 

“j 




( 2 . 1 ) 


where Roman indices refer to the projectile constituents and Greek 
indices to targets. The first two terms in (2.1) are, respectively, 
the kinetic energy and potential energy operators for the projectile 
and are written as 


H - 2 T, + 2 V, . (2.2) 

P j-l J Kj 1J 

Similarly, the third and fourth terms are the target kinetic and 
potential energy operators. The last term is the interaction term 
between ath and jth constituents of the target and projectile. One can 
decouple the center of mass motion of the projectile as 


P 2 
P 


H - 


+ h 


2M A 
n p 


— P — 

with P - 2 P. 

P J-l J 


(2.3) 


where h is the internal Hamiltonian of the projectile which depends 
P 

neither on nor on its canonically conjugate position variable. 


Similar results obtain for the target 


P 2 
T 


h t “ 


“r 


2 P 


(2.4) 


2M n *T ‘ a-1 

where h^ is the internal Hamiltonian of the target. Then the full 

Hamiltonian can be written as 


P 2 ( A + a t ) — „ 

H- + P - - k 2 +h+lu. + V ; ft - 1 (2.5) 

2 M n<VV “nW P 

where the overall center of mass momentum operator is 

F-? p+ 7 x 
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and the projectile momentum compared to the overall center of mass is 


A P 


k - P 


V 


(2.7) 


( 2 . 8 ) 


p ( VV ( w T 

and the interaction is the sum of two body interactions 

\ 

V I ' = V «j 

Ot-l j-1 J 

The first term in (2.5) is the N-body center of mass motion 
energy, decoupled from the other terms. The second is the kinetic 
energy of relative motion of the projectile and the target. The 
projectile relative position variable appears only in the interaction 
term . The projectile and target internal Hamiltonians h^ and h^, 
are coupled to the relative motion through the interaction V^. As the 
separation between the projectile and target becomes larger, tends 
to zero. We assume that well defined states are prepared in the 
entering state and observed in the final state. We define these (in 
operator notation) to be eigenstates of the free projectile- target 
Hamiltonian 

(2.9) 

The full wave function satisfies the Schrodinger equation 

H^ - Erf> (2.10) 

where yfr consists of a superposition of a free state plus a scattered 
state 


(H +H t ) * - E* 
P T 


^ + i> 


scat 


( 2 . 11 ) 


i> - GT * 
scat 


with 

where the Green's function is defined as 


(E'Hp"H<p) *^G-1 


( 2 . 12 ) 
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and the transition operator is defined as 

T - V+VGT (2.13) 

and the wave operator Q as 

- ty (2.14) 

The wave operator satisfies the Lippman- Schwinger equation 
0 - 1+GVfl (2.15) 

so that the transition operator is formally given by 

T - VQ (2.16) 

The goal is to find a series for T. For nucleus -nucleus 
collisions this series was derived by Wilson.^® For a single 
projectile, this series reduces to the familiar Watson-^ series. Using 
the Eikonal approximation, Glauber ^ theory is recovered from 
Wilson's. 

Wilson^ ,40 observed that the transition operator for the 
scattering of a constituent from the jth constituent can be written as 

t.-v.+vGt, (2.17) 

<*j <*j <*j 

which satisfies a Lippman- Schwinger type equation. The wave operator 
which transforms the entering free state up to the collision of a and j 
constituents can be written as 

W 4 ” 1 + 2 G t„ w (2.18) 

“J «kv(«j> <Jk ^ 

The interpretation of the above equation is as follows . The 
propagation to the time just before a and j constituents collide is the 
sum of an operator which brings the initial free state plus the 
scattered part from the scattering of all other fi and k constituents. 
The full wave operator then consists of the wave operator which 
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transforms the system to the a-j collision plus the contribution due to 
scattering of the ath and jth constituents i.e. 

0 - v . + G t . v . (2.19) 

a j “j “j 

which can be written as 

0-1+ 2 Gt.w. (2.20) 

. aj aj 

aj J J 

Wilson then proves that the series given by equations (2.17)- 

(2.20) constitutes an exact representation of the scattering process 

defined by these equations. Consider the product 

V.O-v.w.+v.Gt.w. 
aj aj aj aj aj aj 

-(v +v Gt.)w 
aj aj a j ' aj 


t w 

«j «j 


( 2 . 21 ) 

Summing over a and j one obtains 

( 2 . 22 ) 

“j 

This completes the proof. ^9 

The Green's function G are true N-body operators. One neglects 
binding effects at high energy and replaces G by free N-body operators 


T - 2 V 4 fl «- 2 t. w 

aj “J aj “J “■! 


G^ which satisfy 

(E - 2 T 4 - 2 T ) G - 1 
.1 a o 

J J “ 


(2.23) 


Watson's form of the impulse approximation consists of writing t^. as 


t - v + v G t . (2.24) 

aj aj aj o aj 

so that the above operator acts as a two-body transition amplitude. By 

iteration of the above the multiple scattering series obtains 

T - 2 2 t . + 2 t„, G t . + . . . (2.25) 

a j ‘1 
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which constitutes a formal solution to the exact scattering problem. 

The replacement G-K3 q renders the essentially two-body operators 
and the series for T above, equation (2.25) becomes a series of two- 
body operators. 

The above series reduces to the Watson^l series when the 
projectile is a single particle, as previously mentioned. Next we 
shall derive an optical potential operator whose Born series is 
equivalent to the multiple scattering series expansion (2.25). Such an 
operator is V 0 p t , defined from 


as 


T - V + V G T . 
opt opt opt opt 


V -2 2 t , 

°pt o j «J 


(2.26) 


(2.27) 


From which we obtain 


T - T - 2 t , G t . _ . . . 

opt a j aj aj 


(2.28) 

Retaining the first term in (2.26), the optical model is obtained with 


T _ T “ V G V /(A ‘A—) 
opt opt opt' p T 


(2.29) 


since t 


aj 


opt 

(A p A T ) 


where A ,A^ are the mass numbers of the 


projectile and the target, respectively. 

The approximate Lippman- Schwinger equation for the effective 
potential operator is given by 
0' - 1 + G V 

opt 


where the first order correction to the model is 0 ( ) . 

^ / 

P 


A d A T 


(2.30) 
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II. b Coupled Channel Equations 

Using the multiple scattering theory of nucleus -nucleus collisions 
following Reference 39-42, a set of coupled equations will be derived. 
An equivalent Schrodinger equation will be extracted from the formalism 
developed so far. We shall focus on the Green's function 


G - 2 


k,m,/i 


| g p , m^T , M^k> ^p.m 5 !,/! * k | 
P T 

E - E - E - £• + in 
m n k 


(2.31) 


where E 


P 

m’ 


refer to the projectile and target internal states and 


<VV k 1? 

v — " — 

K 2M A A_ 2m 

n p t 'red 

with the Green's function G from (2.31) inserted into (2.11) and 
projecting onto configuration space yields 


(2.32) 




* “ 8 P.o 


^r.i 


(e T ) 


K <*> 


W 2 r 

(A +A t ) m,M ^ 


,3- ,3- A y? exp(ik |x-y | ) 

d y d d *t 


4w x-y 


9 * 


x ®p,m^p^®T,u^T^®p,m ^p* ^T.u 


x V opt (y,$ p ,? T ) (y,C p .^ T ) 


where 


(2.33) 


V x) “ 


ik»x 


( 2 *) 


3/2 


(2.34a) 
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“ * r i> r 2 '- • • r A } 

P 


(2.34b) 


(s 1( s 2 ,. . 


3 A • 

a t 


(2.34c) 


and 


2 2 
k z - k z + 
m n 


-» VA 
( W 


(<E - ) + (E. - E* ) ] 

p , o m T , o m 


(2.35) 


and gp m (?p) and &p ^(£ T ) are t ^ ie projectile and target internal 
many-body wave funct ions . Now, following Foldy and Walecka^ > we assume 
that the energy transfer is small compared to the incident kinetic 
energy 

k^-k (2.36) 


Using the closure approximation, (2.33) can be written as 


v«*> - «>.. <V •*.. <v >k <*> 

4ir J |X- y) 


2M A 


nVr 


opt V ‘‘T 


(y. i n , $ T > V (y. O f T )... 


p’ 'T' 


(2.37) 


The equivalent Schrodinger equation follows from the above 
equation 

- - 2M A A_ — - 

(V x+ k ) *• <*, < p . < T ) - T75 ^ v opt (X, {p . ? T ) r (». ?p , c T ) 

P (2.38) 

We express the fact that the projectile and target internal wave 

functions are not eigenstates of the optical potential operator and the 

initial states are mixed into various modes of final excited states as 

follows 

V * l % <*> 8p,m <V g T, ( . <V <2 ' 39) 
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The coupled equations then can be written as 


42 


r2 . ,T2, t 2M n A p A T 


(? + k ) ^ (x) - 

x mp 


<VV 


X 

m ' n 


V , , (x)tf , , (x) 

mp.mV m V 


(2.40) 


which is the desired equation. Here 

V ,.'.' 00 - ( e p ,« S r , P l V optl S P.»' 8 T,«' > 


(2.41) 


The coupled equations (2.40) for composite particle scattering 
relate all the entrance channels of the system labelled by projectile 
quantum number m and target quantum number p to all the exit channels. 
In (2.4(1), Ap and are the mass numbers of the projectile and target, 
subscripts m and /j label the eigenstates of the projectile and target; 
M n is the constituent nucleon mass, k is the projectile momentum 
relative to the center of mass, x is the projectile position vector 
relative to the target and 


V (x ’ V v 


2 . t . (x , x.) 
aj aj a’ j' 


(2.42) 


The internal coordinates and have been defined in (2.34) and the 
transition amplitude have been introduced in (2.17). 

Next we shall write the coupled equations (2.40) in matrix form. 
Introducing the wave vector 


(it) 


*00 

*01 <*> 
*10 ( *> 
*11 


(2.43) 
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and the potential matrix 


U (x) 


2M n*p*T 

< A p +A T> 


V 00,00 

v 

v 01,00 

v 

10,00 

v 

11,00 


(X) 

(x) 

(X) 

(X) 


v 00,01 (x) 

v 01,01 (x) 

V 10,01 (x) 

v il,01 (x) 


v 00,10 

v 

01,10 

v 

10,10 

v 

11,10 


(x) . 
(x) . 

(X) . 

<3 • 


( 2 . 54 ) 


The coupled equations (2.40) can be written in matrix form as 




U (5T) * (x) 


(2.45) 


Using the definition of the potential as given by equation (2.42), 
(2.41) can also be written as 

V m/i,mV (x) “ o j <g p.“> ( V g T, ^T* I Vj (X a ,X j ) l g p,m' 


S f P T u . (r ) t . (x ,x ) p ,(r ) d 3 r d 3 r (2.46) 
. J L,pp a aj a j p,mm j j a 




where 


w ( 7> - j *p,» <v $3 <7- { p. 3> «»..• <v d3 «i 

and 


(2.47) 


p T,pp ,( ' r J “ J ®T,m 5 (r a* ^T,a ) ^T )d 


(2.48) 


This completes our review of multiple scattering theory. We shall 
now discuss the formulation of problems of momentum transfer and 
excitation energy deposition in nucleus -nucleus collisions, guided by 
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the multiple scattering approach. Before we proceed, however, a brief 
review of the optical model is in order, since we shall use the optical 
potential in our work. 


II. c The Optical Model 

The optical model is defined as the approximation of (2.40) for 
the elastic scattered part as^’^0 


(v Z + k Z ) * (x) 
where 


_ ^nVr 

( W 

*<*) “ ^00 


W(X) tf(x) 


(x) 


and 

W(x) - (g^ 6p j ol^opt I ®T, 0 «p,0 ) 


(2.49) 


(2.50) 


(2.51) 


with coupling to various excited internal states neglected. This is 
correct at small momentum transfer or near forward scattering. The 
corresponding approximate wave function is called the coherent 
scattered wave and it dominates the forward scattered component. To 
evaluate the optical potential we calculate the Fourier transform of a 
single term of (2.51) 

(&T^o g P,ol t 0** K,q *l 8 T.O ®P.O* ” F T,0 (q) F,0 (q)t y9i (K,q) ( 2 -52) 


where this term will be recognized as the single scattering term of the 

multiple scattering series; q is the momentum transfer and Fp q(<?) and 

— ► 

F t ^(q) are the Fourier transforms of the single particle density of 
the projectile and target, respectively. The nucleon-nucleon 
interaction will be assumed to be constituent averaged as follows^® 
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[N N_ t + Z Z_ t 
p T nn p T pp 


+ (Z L + N Z ) t ] 
p T p T np 


where N , N_ are the neutron numbers and Z , Z_ are the proton numbers 
p T p T r 

the projectile and target, respectively. The optical potential is 
obtained by evaluating (2.42) (summing over constituents) as^® 

W(iT) - I d 3 f p T (f> f d 3 r p p (x + 4 + C) t (k, f) (2.54) 

where pp, p T are the projectile and target single-particle densities 

t(k,y) is the energy and space dependent two-body transition amplitude. 
We have arrived at the expression for the optical potential used to 

analyze heavy- ion scattering cross sections. The input parameters are 

the nuclear ground state densities and the two-body transition 

amplitudes which will be discussed in Chapter IV. 
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CHAPTER III 


MOMENTUM AND ENERGY TRANSFER IN HEAVY- ION COLLISIONS 
In this chapter we shall address the problem of momentum and 
energy transfer in nucleus -nucleus collisions in the framework of 
multiple scattering theory. In the last chapter, we reviewed the 
multiple scattering series, derived the coupled channel equations and 
discussed the optical model. We shall use these tools in the formula- 
tion of the problems of momentum and energy transfer in nucleus -nucleus 
collisions. Starting with the reduced Schrodinger equation for the 
combined system, an approximate expression for the wave operator will 
be derived. Momentum transfer to the projectile as well as its excita- 
tion energy due to collision will then be evaluated. 


III. a Momentum Transfer 

Our starting point will be the formulation of the problem of Ap 
projectile constituents colliding with target constituents as in 
equation (2.1). Following equations (2.1) -(2. 4), the combined 
Hamiltonian has been written in (2.5) as 


H 




<vv 


2 M n<VV 2M nVr 


+ p T k 2 + h + h_ + V, 

D I I 


(2.5) 


The N-body center of mass motion energy, the first term in (2.5) has 
been decoupled from the other terms. The second term is the relative 
motion kinetic energy, with k as in equation (2.7). The third and 
fourth terms, hp and h^ are the internal Hamiltonians of the projectile 
and target with 
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and 


h g (| ) -eg (£ ) 

P P. m P P» m P>® P 

h T g T,M ( ^T ) “ £ T,/i 


(3.1a) 

(3.1b) 


where g^ m (^P) > ^(£^) are t ^ ie i nterna ^ many-body wave functions of 

the projectile and the target respectively; c ^ ^ and c ^ ^ are the 


corresponding eigenvalues. 

The eigenstates 

are orthonormal and com- 

plete 






l 

m 

8 p,»<V 

g P i 


- ‘«T-V 

(3.2a) 


— 





l 

g T#i (C T ) 

g T,; 


- - q) 

(3.2b) 






In 

the overall center 

of mass frame 








p 

IV 

k 



(3.3a) 

p 











P T 

-£ p „- 

-k 



(3.3b) 


a 











and P 

P 

+ p t - 0 




(3.3c) 


Then, the Hamiltonian is 


H - 


1 k 2 + h + h_ + V ; n . 

p T I ’ *red 


n» n A At 

- __E E - reduced mass (3.4) 


2 "red ' (A p + V 

The Schrodinger equation (in coordinate space) can now be written as 


1 

2 H 


*7 2 + h + h_ + V_ 
x p T I 


(X, e p , £ t ) - E v (X, e p , tj) (3.5) 


We seek solutions of the form^^ 
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4> (x,C p ,£ T ) 8 p m (Cp) 


where the bracketed term <f> (x,^ p ,^ T ) J is the desired wave operator. 

Assuming that the source of the beam is along the direction -z, we 
expect 


lim 

Z -+ -oo 


<t> (*.< pI < T >"| - l 


(3.7) 


lim 

Z - 00 


p- — ^ — ► 

|v»(x, Cp. “ 8 p m (?p) ^ (^ T ) 


(3.8) 


Inserting (3.6) into (3.5) yields 


1 -*2 1 ■* i "*2 

_ V - ik*V + ___ k + c + K +e_,+K^ + V 

2 "red X 2 ^red * 2 "red P ° ? ° ^ 


* (x. £ , 4 ) - E * (x,* p ,? T ) 


(3.9) 


where (3.1a), (3.1b) have been used. K p and refer to internal 
motion (Fermi motion) of the projectile and target respectively. We 
now note that 


1 ~*2 

— k Z + e + e_ 

2n P, ° T, ° 


(3.10) 


Solving for the wave operator [ <f> (x,£ p ,£x)] ^- n (3-9) is in gener- 
al difficult. Considerable simplifications occur at high energy when 
[ <f> (x,£ p ,£x)] is assumed to be slowly varying. Then, the second 
derivate in (3.9) can be dropped. ^ Internal Fermi motion K p and K-j 
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can also be neglected. Using (3.10) and the above assumptions, the 
approximate wave operator [ <f> (x.fp,^)] can be written as 


<t> (x. 


i • f > - 


exp 


- P 

k J -® 


2 _ 


I <*'• V V dz ' 


(3.11) 


- exp 


where 


J z 

V 

- an 


(*'. i, £_) dz' 


k - /iv 


-00 I ■»p’ '»X / 


v - relative velocity 


(3.12) 


(3.13) 


Now (3.6) reads 


i> (x, f pt ? T ) 


(2ir) 


3/2 


exp 


- /.I v i <*■. C v dz ' 

IT * 


* 8 P.- ( V s r.p ( «T ) * 


ik»x 


(3.14) 


Using (3.14), we shall next evaluate the total momentum of the 
projectile. Total momentum can be defined as the expectation value of 
the sums of single nucleon momentum operators of the projectile between 
(3.14), i.e. 


P total” <1 * (x ’ ^p’ 


P 

-H T 

j-i Vj 


i> (x, £ p , £ T )> 


< &r 


P 

il v 


,m <V gr M (^)| e xp (-ik.x) exp |+i/v£ V I (x '' V *T )dz ') 

exp |.i/v£ V I^ - V V dZ ) 
exp (ik-x) | g pm ,(? p ) St iM> (^ t )> 


j“l ^P.J 


(3.15) 
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Equation (3.15) is our expression of interest. Note that the 
single nucleon momentum operators act with respect to the internal 
coordinates of the projectile. Note also that we write the interaction 
term Vj(x,^ p , £ x ) as the sum two-body interactions between the ath 
and jth constituent. Since the projectile and target internal 
eigenstates gp m and g XjU are not eigenstates of the wave operator, the 
above is still a coupled channel problem. 

We shall now make a series of approximations in order to solve 
(3.15). Note that the wave operator in (3 . f| ) was obtained from equa- 
tion (3.9) by integration. This assumes that the interaction Vj 
commutes with itself along different points of a straight line trajec- 
tory. If they do not commute, however, then a power series^ solution 
is appropriate as 


4 > (x, ? p , <J T ) 


- 1 v V*’ ♦ ( ;) 2 J. 


Vj (X 1 , $ , € t ) x Vj (X", € , ? T ) + ... 


(3.16) 


- exp 


i 

v 



VjCx', 



e x )dz f 


+ 


(3.17) 


where bracket [ ] + implies an ordering in analogy with time ordering, 

i.e., the argument with earlier z in Vj occurs to the left of the Vj 
with argument with later z. However, if the interactions Vj at dif- 
ferent points along a straightline trajectory commute, then (3.12) 
becomes^ 
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^ (x, £ , | T ) - exp 


(3.18) 


We will assume that the interactions along a straightline trajec- 
tory commute and replace the wave operator (3.17) by (3.18). A plausi- 
bility argument will be given to justify the above. 

When the collision time is short compared to the period of orbital 
motion of the constituent nucleons, then the interactions Vj at two 
successive positions along a trajectory commute, i.e.^,46 


Vj(x', ( e T ) , Vj(x' ' , f p , * T ) - 0 


(3.19) 


Ordering in (3.17) becomes inconsequential, and the upper limit in all 
the integrals can be replaced by infinity to yield (3.18). Using 
(3.18), we can rewrite (3.15) as 

P total- <g p,m ( V g T.M ( Vl exp( ' ik,x) 


Ai/v f v <?, f, ? T )dz'\ -i l 71 
\ P / j-1 f P.J 


^-i/v V ! (x ' ’ V ^T )dz 'J exp (ik*x) | g p,m' ( V T J> 


(3.20) 


In the light of equations (2.43), (2.44) and (2.45) it may be 
instructive to look at the following decomposition of the potential 

matrix^ 


V - V.. + V - V , + V . 

diag off-diag d o.d. 


(3.21) 
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where diagonal and off-diagonal elements have been indicated explic- 
itly. It has been shown that for forward scattering at high 
energies, the diagonal elements in V dominate over the off-diagonal 
elements^® 


V.. » V __ .. 

diag off-diag 


(3.22) 


Then (3.20) can be written as 
P tot.r <s p,m <y s T,p«T ) l exp flk-x) exp ( I/v J 


(V d + v „.d> dz 


') 

(' iA l.^ + Vd )d2 ) exp<+lk -*>|y»’<V *!,*« )P 

I J 

- < gp i m ( ^p ) g Tl M ( ^T ) l exp ( i/v £ V d ( ^’ V V dZ ) (1 + ;£Vd dz ' + 


A 

P -*» 

-il V 

j-1 *P,J 


P 

-iZ 7 
j-i 5 p.j 


exp 


(• I/v J.?d (? ’ v v dz ’) (1 -;J. 


00 _ 


V . dz'+. . .) 
o.d 


X | g p,m'^p^ 

■ %,» <V jjV*'- V V dz '^ 

exp (-1/vJ V, r T )dz.) | g 


(3.23) 

A 

P 

■iz r 

j-i ^p.j 

(3.24) 


+ other terms . 


It is now obvious that the other terms in the series will involve 
the off-diagonal elements which are smaller (by (3.22)). The first 
term in the series involves only diagonal elements. In the spirit of 


/ 
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the coherent approximation^ 2 , we shall retain only the first term in 

( 3 . 24 ). 

Now, applying the identity^ 

2 

2 


e iA B e' iA - B + i [A, B] + L [A, [A, B] ] + . . . . 


( 3 . 25 ) 


leads to 


P total“ <g p,m ( V 8 T,m ( ^T ) 


A 

-i l V 


j-1 S P-j 


g p,m ( V 


P 


<g p,m <V g T, / i (C T ) E J>~ ’ V " T 


j-1 P.J 


V„(x', *t> — 

V 


®p,m ( V St./V* 


( 3 . 26 ) 


Notice that the first term can be thought of as the momentum of the 
projectile before the collision. It may be taken to the left-hand side 
of (3.26) and momentum transfer defined as 

p p _ p 

total before transfer 


<g p,m ( V g T,M (? T ) 


j"^V d (x' , £ p . £ t ) 


j-1 P.j 


dz* 

v 


dz' 


*p." «p> g T.i. << T )> 

^ °° -*■ -» 

<g (? ) &r (£—) I 7 i f V nn f ( - X ' • O ^ — - 

6 p.o v V 6 T,o vs T / C p j J.* °P C P T v 

g p,o ( ^p ) 


( 3 . 27 ) 


( 3 . 28 ) 


where in the spirit of the optical model discussed in the last chapter, 
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* ^ oq has been used in going from (3.27) to (3.28). 


Having identified through equation (2.27) the optical potential 


opt 


A™ A 

r x p 

o-l j-1 




(3.29) 


we can substitute the above in equation (3.28) and obtain 

A ^ A_ A 

-»■ -*■ -*• | P ■* I 

^transfer ” < ®p , o ^p^ ®T,o^T^ 


r r i i 

j-i S.j J -« «-i j-i 


'.J <*■■«,. v * 1 


h.o <V g T,o (? T )> < 3 ' 30 > 


Equation (3.30) above is the desired expression for momentum 
transfer in the relativistic heavy ion collisions. Using the multiple 
scattering theory formulation of nucleus nucleus collisions of 
Wilson, 39,40 we have seen that the optical potential can be calculated 
in terms of transition amplitudes rather than two-body potentials. We 
obtained above an expression for momentum transfer in relativistic 
heavy ion collisions using these experimentally well determined transi- 
tion amplitudes. Equation (3.30) has been explicitly evaluated in 

the Appendix, with the result expressed as and P^_ - 

r tran tran 

i i 

£ . where P^ J is the momentum transfer in the collision 

aj tran tran 

between nth and jth constituents. Explicitly, these are given by 


KLn -/ d % VV f W 


— *• f <*> a , 

V, I t . (x +£-£„) — 

£P J -» P v 


(3.31) 
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where a and j refer to the ath and jth constituents of the target and 
projectile, respectively; and Pj(£p) are their single particle 

densities, and the gradient is with respect to projectile coordinates 
only. The total momentum transfer is 

P - - 2 2 

tran , tran 

® j 


- - <Vt> I ‘Vy V j d3 «T ' T «T> 


n 


(x', ,e T ) 


dz ' 


v 

(3.32) 


where p (£p) and p are the nuclear densities of the projectile and 

target with mass numbers Ap and respectively, and the two -body 
transition amplitude t . has been averaged over constituent type as 


t - (ApA^ 

with 2p, Z’l the projectile and target charge numbers; Np, Nx the 
neutron numbers and Ap, Ax the mass numbers, respectively. (The 
notation for the transition amplitudes is tpp for proton-proton, t nn 
for neutron -neutron, and t n p for neutron-proton.) 

A unique feature of this work is the use of the two-body transi- 
tion amplitude (in (3.32)) which is complex. The constituent averaged 
two-body transition amplitude is obtained from the first order t-matrix 
used in previous studies of nucleus -nucleus collisions as 


Z L t + N N t + 
p T pp p T nn 


(N Z + N Z ) t 
p T T p np 


(3.33) 
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•3/2 


-*2 

t (e, x) - - +1 e/m a(e) [a(e) + 71 r2wB(e7j exp(- ) 

* ‘ ‘ L— ' An / \ 


2B(e) 


(3.34) 


where e is the two-nucleon kinetic energy in their center of mass 
frame; a(e) is the nucleon-nucleon total cross section; a(e) is the 
ratio of the real to the imaginary part of the forward scattering 
amplitude, and B(e) is the nucleon-nucleon slope parameter. Values for 
these parameters, taken from various compilations, are listed in refer- 
ences 39-42. 

The resulting momentum transfer, equation (3.32), where the tran- 
sition amplitude (3.34) appears is also complex. The real part of the 
momentum transfer, which comes from the real part of the complex opti- 
cal potential, is the contribution arising from elastic scattering. It 
is purely transverse. The imaginary component, which comes from the 
absorptive part of the complex optical potential, is the longitudinal 
momentum transfer. To show this, we symbolically rewrite equation 
(3.32) as 


'transfer " <«R + 1 V b 


(3.35) 


where i «\/-l and b is the unit vector transverse to the beam 

A A 

direction. If z denotes the beam direction (with -z the source) then 


we know that 

ib - - z (3.36) 

since i -V-l is an operator which rotates a unit vector counter 
clockwise through n/2 radians . ^ Therefore (3.35) becomes 
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8 transfer ' \ b ' Q I 2 < 3 - 37 > 

which can be relabelled as 

transfer ’ ^ > ' »| * < 3 - 38 > 

The transverse (Qj_) and longitudinal (Q|) components can be immediately 
written as 


^ • Vr . 

3 -* 

d ? p (e ) 
p p p J 

d 3 ? x P T (^> 

V 

f 00 _ ^ d z 

Re t (x' + i - $ ) _ 

* -co ** V 

and 




(3.39) 

V - Vt . 

d u p a ) 
p p p . 

d 3 £ T p t (C t ) 

\ 

f ® -*• A 7 

[ Im t (X' + S - $ ) 

-oo ^ v 


(3.40) 


The above has a close analogy in optics where complex refractive 
index is used to account for absorption. The resulting attenuation of 
the incident wave is then along the direction of incidence. Similarly, 
momentum loss due to inelastic scatterings and absorption as the pro- 
jectile traverses through the target is in the longitudinal direction, 
while the elastic scattering gives rise to transverse momentum transfer 
(specially true at high energy where forward scattering dominates). 

Illb. Energy Transfer 

Having developed the formalism for momentum transfer, we shall 
also find an expression for energy transfer (excitation energy) to the 
projectile. Again, our basis will be the multiple scattering theory 
framework discussed previously. We will assume in light of the discus- 
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sion in Appendix B that the sudden approximation is valid, i.e., colli- 
sion time t co n is short compared with the period of orbital motion of 
nucleons in the nucleus. This is best evaluated by taking the expecta- 
tion value of the sum of nucleon energy operators (the internal 
Hamiltonian of the projectile H p ) with approximate wave functions given 
by (3.14). 

A 


E - <rf> (x,f £ )| _ _L 2 ~V 2 + 2 V.. 

P T I 2m n j-1 Vj i<j 1J 


(x,C p ,^ T )> (3.41) 


Expanding the above yields 


<S p,m «p> S T, M <€ T ) l eXp <* lk * x > exp^i/v J Jj(x» , | p , ^ )dz j 

exp ^-i/v | V r (x', 7 p> 


1 ^ -*•2 
_ 2 Vt + 2 V 

2m n j-1 S.j i<j J 


exp (Ik-x) | g p „,(S p ) 8 t , p .(« i ) > 


‘ <s p.» <V *t,„ w t>| 


( V V o.d )d2 


■) 


2 V 


2m j-1 
n J 


*P.j 


+ S V. . 

i<j lJ 


exp 


(-W | 


(V d +V o.d )d2 ' 


) g p,m' ( V ®T.M 


<g p i o ( V g T, 


exp 


( i/v SI v 


(X', ? T )dz')x 


*p’ 
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The explicit derivation of the above expressions is detailed in 
Appendix A. The final result for energy transfer is 

E* - E - E +0 (Jl) 
o A 

P 




Evaluation of the above equation for E* yields an expression for energy 


transfer as 


* i f 3. - r 3- 

E - — - ApA T J d £ p P p (£p) J d 

n 

(3.45) 



where the definition of V 0 p t , equation (3.29) has been used, and we 
neglect terms 0(1/Ap). All the symbols above have been defined in our 
previous derivation in section Ilia. 


CHAPTER IV 


NUMERICAL RESULTS ON MOMENTUM AND ENERGY TRANSFER DISTRIBUTIONS 
In this chapter, we address specifically the momentum and energy 
transfer distributions in nuclear fragmentation reactions. Having 
studied in Chapter II the multiple scattering theory of nucleus nucleus 
collisions and in Chapter III the formulation of the problem of energy 
and momentum transfer in heavy ion collisions based on the multiple 
scattering theory, we shall discuss a variety of topics in heavy ion 
reactions with these insights. First, we briefly review the experimen- 
tal situation as it pertains to our topics of interest. We next dis- 
cuss the currently available theories of heavy ion fragmentation 
reactions. We present calculations on the momentum and energy transfer 
distributions in relativistic heavy ion collisions and indicate how to 
incorporate our results in the existing theoretical framework. Our 
primary focus will be on the momentum spectra of fragments, including 
the effects of momentum transfer on the spectra. Input parameters are 
briefly discussed next. Numerical results on momentum "downshifts" and 
widths are presented for a variety of projectile- target combinations 
and compared with experiment. 

IV. a The Experimental and Theoretical Background 
Most of our information on fragmentation reactions comes from 
single particle inclusive reactions of the type^®'^l 
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Beam + Target -+ Fragment + Other Products (4.1) 

Two concepts stand out in describing the projectile (or target) frag- 
mentation data. They are called "limiting fragmentation" and "factor- 
ization," respectively. "Limiting fragmentation" implies that 
distribution of fragments approaches a limiting form as the bombarding 
energy is increased. "Factorization" means that the cross section for 
production of a particular fragment (written as apg'j') is a product of a 

F 

factor 7 x which depends on the target and 75 which depends on the beam 
and fragment , i . e . 


a FBT ” 7 T 7 B 

Further analysis-^ suggests that another prescription, 
factorization" could better describe the data 


(4.2) 

called "weak 


F T 
a FBT “ 7 B 7 B 

T 

where 75 depends both on the beam and the target. 


(4.3) 


Experimental measurements^ Q f widths and means of longitudinal 
and transverse momentum spectra of fragments of and 16 q on targets 
ranging from H through Pb revealed the following features. In the 
projectile rest frame, the longitudinal momentum distributions for all 
fragments from an d 16 q ( w ith the exception of protons) show, irre- 
spective of beam energy and target, the following characteristics: 

(a) Gaussian shape, with root mean square widths o(P||) “ 50 to 
200 MeV/c and downshifts CP y> - - 20 to - 130 MeV/c. Thus the mean 
velocities of the fragments are less than that of the beam. 

(b) Rms widths o(P||) an< * are ®<l ua l to an accuracy of 10%, 

consistent with isotropic production of fragments. 
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(c) a(Pj|) and <P|> are independent of target mass and beam energy 
but depend on the beam and fragment. 

The general trend of the widths a(P||) reproduced by a parabolic 
dependence on the fragment mass where 


a(P||) - 2o q [x(l-x)] l/2 



Fragment mass number 
Projectile mass number 


(4.4) 


where o Q is a constant. 

The parabolic dependence has been explained by a variety of theo- 
retical approaches^" 22 including the conservation of momentum. 

Feshbach and Huang^O and later Goldhaber^ pointed out that the 
Gaussian momentum distributions could be understood by treating the 
fragmenting nucleus as a Fermi gas and assuming (a) momentum conserva- 
tion (b) no correlation among nucleons and (c) neglect of anti- 
symmetrization of the single particle distributions. Based on this 
model, the parabolic dependence of the widths a(P||) on fragment mass Ap 
could be reproduced (as in equation (4.4)), as well as the isotropy 
cr (P |j ) =* a(P^) could be understood. More importantly, Goldhaber^ 
also pointed out that the widths and means of momentum distributions 
could be modified according to 

a' 2 _ ct 2 + q 2 (4.5a) 

x x A 2 

P' - P + £ 0 (4.5b) 

X X . 


where a x is the unmodified Gaussian width, a is the modified width, 
Q x is the momentum transfer in any direction x; F and A are the frag- 
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ment and parent mass numbers respectively and unmodified widths are 
related to Fermi momentum Pp of the parent as 


2 _ F(A-F) 2 


(4.6a) 


F(A-F) P F 
A 5 


(4.6b) 


The longitudinal momentum downshift is given by (4.5b) as 



(4.7) 


where Q|| is the magnitude of the longitudinal momentum transfer 
obtained from equation (3.40). Recalling that Q| is a function of 
impact parameter, an appropriate method for choosing it for each frag- 
ment is necessary. Before that however, we shall present input parame- 
ters required in our calculation. (Note the comprehensive list of 
references 53-88 on the experimental and theoretical overview of the 
subject) . 


IV. c. Input Parameters 

Having derived expressions for momentum and energy transfer, we 
shall use them to calculate quantities of interest. 

The NN transition amplitude in coordinate space has been intro- 
duced in equation (3.34) and the notation explained. 

Values of a(e), a(e) and B(e) are discussed in references 62 and 
86, including an estimate of experimental uncertainty. Average values 
are given by 



1 
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a(e) - 


[<Zp + Z T ) crp p (e) + (N p + N T )<7 np (e)] 


(Ap + Aj) 

— l(Zp + Zf) <7pp(e) Qpp(e) + (Np + N'j) <*np(e) a np ( e ) ] 


(4.8a) 


(4.8b) 


(A p + A T ) (<Tp p (e) + a np (e)] 

Some representative values are listed in Table 1. In addition, 
the "non-diffractive" slope parameters B(e) from these references will 
be replaced by "diffractive" slope parameters 


B(e) - [ 10 + .5 In (S'/S 0 )] (GeV/c) 


-2 


(4.9) 


with S Q * 1 (GeV/c) "2 and S' - square of the NN center of mass ener- 
gy. This is because diffractive slope parameters are appropriate for 
scattering near forward directions. The single particle densities are 
related to the nuclear wave function by 

p k < 7 > - - r <6„<r> I S <r - 7) | g 0 <r» (4.10a) 

p A p a-1 

p ( T ) - J: £ T <g(L) I s (7 - S ) I g (£ t )> (4.10b) 

*T A T j-1 ° 1 J > o i 

for the projectile and target, respectively. The above are understood 
as 


P A ( r > 
P 




'n ( 


r ) + Z 



(4.11a) 


*T 


( s ) - -1 [n t P n ( s ) + z T p ( 7)] 


A-p 


(4.11b) 


where Zp, are the proj'ectile proton and neutron numbers and Zj, N^ 
are the target proton and neutron numbers, respectively. We assume for 
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Table 1. Values of NN input parame ters o(e) for NN - np and NN-pp 
for various incident energies. Data from compilations in references 
62 , 86 . 

E (A GeV) ^ ( e ) (mb) a (e) (mb) 
inc np pp 


.025 

378.255 

130.31 

.050 

140.38 

41.00 

.075 

98.35 

39.63 

.100 

71.73 

27.23 

.125 

59.15 

23.83 

.150 

52.09 

24.13 

.175 

46.59 

23.79 

.200 

42.80 

23.09 

.225 

40.09 

22.45 

.250 

37.94 

22.21 

.275 

36.29 

22.37 

.300 

35.15 

22.78 

.400 

34.03 

25.19 

.500 

34.82 

32.46 

.600 

36.03 

39.28 

.700 

36.98 

43.02 

.800 

37.84 

45.51 

.900 

38.76 

47.01 

1.000 

39.68 

47.65 

2.000 

42.96 

45.18 


3.000 


43.19 


42.50 



light nuclei p n (r) “ pp (r) since Coulomb repulsion by protons plays 
only a small role in light nuclei. Proton densities are taken from 
compilation of charge radii from electron scattering experiments. The 
effect of finite proton radius is taken into account by extracting the 
matter densities according to reference 62. The nuclear charge density 
is taken as 

P CfA ( r) - J P p (r') p a (r + r') d 3 r' (4.12) 

where pp(r') is the proton charge density. 

In this work, the densities of nuclei for A > 20 was taken to be 
of the Woods -Saxon type 



(4.13) 


with R the half density radius and a related to the skin thickness 
as t - 4.4a. For A < 20, Harmonic -well densities were used. The 
parameters are listed in Table 2. 


IV. c. Numerical Results 

A. Momentum Downshifts 

Experimental data on momentum downshifts are available^ for the 
fragmentation of ^0 an< j on targets ranging from H through Pb. In 
Figures 1-18, momentum transfers to ^ 0 ( 2.1 a. GeV) , l^c(2.1 A GeV) and 
(1.05 A GeV) projectiles are plotted as a function of impact param- 
eter b (fm). Momentum transfer is in units of MeV/c. Longitudinal and 
transverse momentum transfers are shown separately. 
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Table 2: Nuclear charge distribution parameters 

from electron scattering data 
(HW - Harmonic -well; WS - Woods -Saxon) 


Nucleus 

Distribution 

7 (HW) or 
t, fm (WS) 

a, fm (HW) or 
R, fm (WS) 

*Be 

HW 

.611 

1.791 

12 C 

HW 

1.247 

1.649 

16 0 

HW 

1 . 544 

1.833 

27 A 1 

WS 

2.501 

3.05 

64 Cu 

WS 

2.504 

4.20 

108 Ag 

WS 

2.354 

5.139 

^La 

WS 

2.354 

5.71 

208 Pb 

WS 

2.416 

6.624 

Harmonic -well: 

p c ( r ) - Po A + y 



1 + exp ft 2 ) 


Woods-Saxon: p c ( r) - p 0 


-1 
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4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 


Impact parameter, fm 

llgtirc I: Momentum transfer (McV/c) to ,r, 0 projectile by the *Be target as a 

function of impact parameter (fm) in the reaction 16 0 (2.1 AGeV) 9 0e -» 
Projectile Fragment + X, tvhere X is unidentified. 1 Iarmonlc*well 
densities were used for both the projectile and tlu» larpei 
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l igurc 2: Momentum transfer (MoV/c) to ^'O projectile by the ,J C target as a 

function of Impact parameter (fm) In the reaction u O (2.1 AGeV) + ,J C -♦ 
rrojeclile Fragmonl + X, where X Is unldenllfled. Harmonic-well 
densities were used for both the prc'ectlle and the target. 
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Impact parameter, fm 

Figure 3: Momentum transfer (MeV/c) to *‘0 projectile by the A1 target as a 

function of impact parameter (fm) in the reaction ! *0 (2.1 AGeV) + AJ -» 
Projectile Fragment -*• X, where X is unidentified. Harmonic* well density 
was used for the projectile and Woods-Saxon density for the target. 
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Impact parameter, fm 

Figure 4: Momentum transfer (MeV/c) to ,8 0 projectile by the Cu target as a 

function of Impact parameter (fm) In the reaction **0 (2.1 AGeV) + Cu -♦ 
Projectile Fragment + X, where X Is unidentified. Harmonic-well density 
was used for the projectile and Woods-Saxon density for the target. 




8.0 8.2 8.4 8.6 8.8 9.0 9.2 


I 

Impact parameter, fm 

Tig. 5 Momentum transfer (MeV/c) to ’*0 projectile by the Ag target as a 

■ function of impact parameter (fm) In the reaction (2.1 AGeV) + Ag -♦ 

Projectile Fragment 4 X, where X Is unidentified. Harmonic-well density 
was used for the pro|ectlle and Woods-Saxon density for the target. 
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Impact parameter, fm 


Figure 6: Momentum transfer (MeV/c) to '*0 projectile by the Pb target as a 

function of Impact parameter (fm) in the reaction **0 (11 AGeV) ♦ Pb -* 
Projectile Fragment + X. where X Is unidentified. Harmonic-well density 
was used for the pro|ectile and Woods-Saxon density for the target. 
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Impact parameter, fm 


Figure 7: Momentum transfer (MeV/c) to l, C projectile by the De target as a 

function of Impact parameter (fm) In the reaction ,J C (1.05 AGeV) + ’Be 
-* Projectile Fragment X, where X Is unidentified. Hatmonlc-well 
densities were used for both the projectile and the target. 
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Impact parameter, fm 

Figure 8: Momentum transfer (MeV/c) to IJ C projectile by the W C target as a 

function of impact parameter (fm) In the reaction tJ C (1.05 AGeV) 
,J C-* Projectile Fragment X, where X Is unidentified. Harmonic-well 
densities were used for both the projectile and the target. 
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5.8 


6.0 


6.2 


6.4 


Impact parameter, fm 

Momentum transfer (MeV/c) to projectile by the J7 AI target as a 
function of Impact parameter (fm) In the reaction ll C (1.05 AGeV) + 

37 A1 -♦ Projectile Fragment + X, where X Is unidentified. Harmonic- 
well density was used for the projectile, and Woods-Saxon was used for 
the target. 
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Figure 10: Momentum transfer (MeV/c) to projectile by the Cu target as a 

function of impact parameter (fm) In the reaction t2 C (1.05 AGeV) + Cu 
-* Projectile Fragment 4- X, where X la unidentified. Harmonic-well 
density was tised for the projectile, and Wooda-Saxon was used for the 
target. 
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Impact parameter, fm 

igure 11: Momentum transfer (MeV/c) to projectile by the Ag target as a 

function of Impact parameter (fm) In the reaction (1.05 AGeV) + Ag 
-» Trojectlle Fragment X, where X It unidentified. Harmonic-well 
density was used for the projectile, and Woods-Saxon was used for the 
target. 
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Fig. 12 Momentum transfer (MeV/c) to n C projectile by the Pb target as a 

function of Impact parameter (fm) In the reaction ,2 C (1.05 AGeV) Pb 
Projectile Fragment + X, where X Is unidentified. Harmonlc*well 
density was used for the projectile, and Woods-Saxon was used for the 
target. 
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4.6 4.8 5.0 5.2 5.4 

Impact parameter, fm 


Figure 13: Momentum transfer (MeV/c) to ,3 C projectile by the De target as a { 

function of Impact parameter (fm) In the reaction (2.1 AGeV) *Be 
-» Projectile Fragment • X, where X Is unidentified. Harmonic-well 
densities were used for both the projectile and the target. 
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igure 14: Momentum transfer (MeV/c) to ,J C projectile by the ,J C target as a 

function of Impact parameter (fm) In the reaction ,J C (2.1 AGeV) + 
n C-* Projectile Fragment X, where X Is unidentified. Harmonic-well 
densities were used for both the projectile and the target. 
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Impact parameter, fm 

Fig. 15 Momentum transfer (MeV/c) to 1J C projectile by the target as a 

function of Impact parameter (fm) In the reaction ,J C (2.1 AGeV) + J7 A1 
-♦ Trojectlle Fragment + X, where X Is unidentified. Harmonic-well 
density was used for the projectile, and Woods-Saxon was used for the 
target. 
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Figure 16: 





7.2 


7.4 


Impact parameter, fm 


Momentum transfer (MeV/c) to ,l C projectile by the Cu target as a 
function of Impact parameter (fm) In the reaction ,2 C (2.1 AGeV) + Cu 
-4 Projectile Fragment + X, where X Is unidentified. Harmonic-well 
density was used for the projectile, and Woods-Saxon was used for the 
target. 
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Impact parameter, fm 

Fig. 17 Momentum transfer (MeV/c) to ,J C projectile by the Ag target as a 

function of impact parameter (fm) In the reaction ,J C (2.1 AGeV) + Ag 
-4 Projectile Fragment + X, where X It unidentified. Harmonic-well 
density was used for the projectile, and Woods*Saxon was used for the 
target. 
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Fig. 18 Momentum transfer (MeV/c) to l, C projectile by the Pb target as a 

function of impact parameter (fm) in the reaction ,J C (2.1 AGeV) ♦ Pb -> 
Projectile Fragment + X, where X is unidentified. Harmonic-well 
density was used for the projectile, and Woods-Saxon was used for the 
target. 
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From the figures, two features are readily apparent. First the 
longitudinal momentum transfer is larger than the transverse, indi- 
cating the primarily absorptive nature of nuclear collisions at this 
energy. Second, the predicted momentum transfers decrease rapidly with 
increasing impact parameter. 

Another feature that stands out is that momentum transfers are 
only slightly target dependent. The choice of the impact parameters 
was of consequence in this and will be discussed shortly. The longitu- 
dinal momentum transfer to ^0 rises from -230 MeV/c for ^Be target to 
-320 MeV/c for 208p^ target at the closest impact parameters (see 
figures 1 and 6) . The transverse momentum transfer increase is only 
-90 MeV/c to -120 MeV/c. The projectile dependence is more noticeable, 
however. The longitudinal momentum transfer to ^ C (at 1.05 A GeV) is 
-155 MeV/c (figure 7) for ^Be and rises to -180 MeV/c (figure 12) for 
208pb. Transverse momentum transfer is -25 MeV/c - 30 MeV/c only for 
all targets. At incident energy of 2.1 A GeV, however, the energy 
dependence of transverse momentum transfer is noticeable. It ranges 
from -60 MeV/c for the ^Be target (figure 13) to -70 MeV/c for 208pb 
(figure 18). Longitudinal momentum transfer hardly changes between 
1.05 A GeV and 2.1 A GeV. The increase in transverse momentum transfer 
is related to a(e) which is the energy -dependent ratib of real to the 
imaginary part of the NN forward scattering amplitude. These trends 
show that while the longitudinal momentum transfer remains fairly 
constant at high energies, the transverse momentum transfer is highly 
energy - dependent . 
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In Figures 19-21 and Tables 3-5, our results on momentum down- 


shifts are reported. Experimental data have been averaged over 
isotopes using 


<Pn> 


ave 

exp 


"FI <P || > F1 + a F2 <P ||> 2 + --- 


°n + a F2 + --- 


(4.14) 


where ap^, <7p2 are the fragmentation cross sections for Isotopes 1 and 
2, respectively, and <P||>F1* <P || > F2 are t^eir corresponding downshifts 
(experimentally observed) . These calculations are presented in Figures 
19-21. To translate the calculated longitudinal momentum 
transfers into "Momentum downshifts," we follow the following prescrip- 
tions . 

(1) The semi-empirical fragmentation code NUCFRAG®^ developed at 
NASA Langley is used to assign a range of impact parameters - b2, 
b2 - bj, for each fragment Ap - 1,2,3. .. (Ap-1) where Ap, Ap are the 
fragment and parent mass numbers, respectively. This range is divided 
into - 30 intervals and corresponding longitudinal momentum transfer 
calculated for these values of impact parameters. An arithmetic aver- 
age is done and the average is multiplied by the Goldhaber factor. 

This is the "momentum downshift" for fragment Ap. Standard deviation 
of the mean is computed following standard procedures. 

(2) An impact parameter b is uniquely assigned by NUCFRAG for 
each fragment Ap where Ap - 1 , 2 . . . (A p -1) . The corresponding longitudi- 
nal momentum transfer is calculated. The average longitudinal momentum 
transfer is taken as the same as the above. This average, multiplied 
by the Goldhaber factor (Eq. (4.7)) is the "momentum downshift" for 
that fragment. The results are compared from the two approaches. 
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Fragment mass number 

Figure 19: Momentum downshifts of projectile fragments in the reaction * } C 

(1.05 AGeV) ■*• Target -4 Projectile Pragment 4 X, where X Is 
unidentified and targets are *Be, ,J C J 'A1, Cu, Ag and Pb. Harmonic- 
well densities were used for A < 20, Woods-Saxon for A > 20. 
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Fragment mass number 

Fig. 20 Momentum downshifts of projectile fragments In the reaction ,J C 

(2.1 AGeV) Target -* Projectile Fragment * X, where X is unidentified 
and targets are ’Be, ,2 C, 27 AI, Cu, Ag and Pb. Harmonic-well densities 
were used for A < 20, Woods-Saxon for A > 20. 
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Figure 21: Momentum downshift* of projectile fragments In the reaction «0 

(2.1 AGeV) ♦ Target -> Projectile Fragment + X, where X Is unidentified 
and targets are *De, IJ C, 17 A\, Cu, Ag and Pb. Harmonic-well densities 
were used for A < 20, Woods-Saxon for A > 20. 
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Extensive numerical work revealed that the second approach, while being 
more efficient yields results comparable to the first approach, not 
differing by more than 7 MeV/c. When accuracy was desired, the first 
approach was always preferred over the second. 

It may now be worthwhile to discuss the experimental results and 
theoretical calculations. It is clear from figures 19-21 that the 
momentum downshifts are overpredicted in all cases except for a few. 

The source of this discrepancy is the impact parameters obtained from 
the geometrical, semi -empirical code NUCFRAG. The projectile and 
target densities are approximated in the above code as uniform spheres 

with R - 1.26 aV 3 (fm ). This is obviously an oversimplifi- 
cation. Electron scattering®^ from nuclei reveal that nuclei possess 

diffuse surfaces. Realistic charge distributions (Woods -Saxon for 
example) take the diffuseness into account through skin thickness 
t(fm) , which is a measure of the distance where nuclear density falls 
from 90% to 10% of its value. A sharp cut-off of the density thus 
neglects the extended, diffuse nuclear surface. Since realistic nucle- 
ar densities were used in the calculations of longitudinal and trans- 
verse momentum transfer, it is likely that the impact parameters from 
NUCFRAG are not very realistic. An alternate procedure would be to 
make our calculations compatible with uniform density calculations from 
NUCFRAG. 

This was accomplished in two steps. Realistic densities were 
replaced by uniform densities £ R- 1.26 A^/^]in the momentum transfer 
calculations. The zero-range of the two-body interaction in NUCFRAG 
was implemented in our calculation by reducing the range arbitrarily by 
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^ / g \ 

a factor of five i.e. B'(e) - _ where B(e) was given previously by 

5 

equation (4.9) all proj ectile- target combinations. This amounted to 
B'(e) -.08 fm"2 at 2.1 A GeV, for example. The effect on our calcula- 
tion can be seen in figures 22*24. In figure 22, the imaginary part of 
the optical potential V 0 p t (r) for . 9g e collision at 2.1 A GeV is 
plotted using Harmonic-well densities as well as uniform densities with 
variable slope parameter B(e) . The first-order optical potential 
follows closely the density distribution of the nuclei (actually ex- 
tends beyond due to finite range of the interaction) . As uniform 
density was substituted and the range of the interaction reduced by a 
factor of 5, the shape of the potential became steeper, thus approxi- 
mating a sphere with a sharp cutoff. The effect on the longitudinal 
momentum transfer to ^0 i n the collision of the l^O-^Be pair (Ei nc - 
2.1 A GeV) can be seen in figure 23. Compared with the realistic 
density calculations (see also figure 1), the magnitude of the momentum 
transfer is reduced. Setting B(e) - 0 outright involved numerical 
difficulties . 

The momentum downshifts of Oxygen fragments are plotted and 
compared with experimental data^ in figure 24. The experimental data 
have been averaged according to (4.14). Notice the significant agree- 
ment between theory and experiment. Similar results obtain for 
fragmentation at 2.1 A GeV and 1.05 A GeV. 

It may be appropriate to point out that the above procedure can be 
reversed and theoretical calculations utilized as an impact parameter 
"gauge". With experimentally observed downshifts as inputs, one can 
calculate with realistic densities the impact parameters where the same 


Imaginary Vopt (r) f MeV 
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l'ig. 22 Imaginary part of the optical potential Vopt(r) In MeV as a function of 
separation r In fm for Oxygen with Beryllium at 2.1 A GeV. Harmonic 
well densities were used for both the projectile and the target and 
compared with uniform density calculations. Slope parameter was 
modified to B'(e)-D(e)/2 and B*(e)»B(e)/5 respectively to simulate a 
zero-range interaction. 
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4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 


Impact parameter, fm 

rig. 23 Longitudinal momentum transfer In MeV/c to Oxygen by the 

Beryllium target as a function of Impact parameter In fm. Harmonic 
welt densities were used for both the projectile and the target and 
compared with uniform density calculations. Slope parameter was 
fnodlfied to D*(e)-B(e)/2 and B*(e)-B(e)/5 respectively to simulate a 
zero-range Interaction. 
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Fragment mass number 

Fig. 24 Momentum downshifts of projectile fragments In the reectlon '*0 

(2.1 AGeV) 4 - Target -*■ Projectile Fragment + X, where X Is unidentified 
and targets are ’Be, lJ C 27 Al, Cu, Ag and Pb. Uniform densities were 
used for the projectile and the targets. Slope parameter was modified 
(B'(e) - B(e)/5) to approximate a zero* range Interaction. 
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downshifts obtain. The practical advantage of these impact parameters 
is that they can be utilized in localizing a reaction event as well as 
used as guides to other calculations such as those of transverse 
momentum transfer, fragmentation cross-sections, and Coulomb dissocia- 
tion cross - sections . 

The energy dependence of the momentum transfer is a quantity of 
fundamental interest. As noted previously, the NN input parameters are 
energy dependent. Longitudinal and transverse momentum transfers have 
been plotted for - ^Be pair at incident energies Ej nc <-0.2, 0.4, 
0.6, 0.8, 1.05 and 2.1 A GeV in figures 25 and 26, respectively. Note 
the gradual increase of longitudinal momentum transfer as the bom- 
barding energy increases, attaining a "limiting" value at E^ nc > 1 
A GeV. The transverse momentum transfer can be readily obtained from 
the longitudinal momentum transfer by multiplying the latter by a(e), 
where a(e) is the ratio of the real to the imaginary part of the NN 
forward scattering amplitude. As is well-known, a(e) < 0 at 2.1 A GeV 
which implies a repulsive, real part of the optical potential. The 
corresponding transverse momentum transfer is from a repulsive mean 
field and is positive (according to our sign convention). At lower 
bombarding energies, V 0 p t (r) - -(V^ ea l + iVj ma g) holds and the corre- 
sponding mean field is attractive. Transverse momentum transfer is 
negative indicating that deflection to negative scattering angles is 
feasible. Coulomb effects have been ignored in our calculations. The 
above predictions should be verifiable in sophisticated experiments in 


the future. 
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Impact parameter, fm 

Figure 25: Longitudinal momentum transfer (MeV/c) to projectile by the y Be 

target as a function of Impact parameter (fm) In the reaction 
1s CXE (ne )4- 9 Be — >Pro)ectile Fragment-fX. Incident energies 
E, ne -2?l,t .05,.8,.6,.4 and .2 A GeV were used. 
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f igure 26: Transverse momentum transfer (MeV/cj to projectile by the *tJe 

target as a function of Impact parameter (fm) In the reaction '*0 (B| na ) + 
’Be -♦ Projectile Fragment 4 X. Incident energies E |M - 2.1, 1.05, .8, .6, .4. 
.2 AGeV were used. 
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B. Momentum Widths 

Equation (4.5a) is Goldhaber's result for the modification of the 
fragment momentum width due to collisional momentum transfer. Specifi- 
cally, the following equations can be derived from (4.5a) 



(4.15a) 


(4.15b) 


where a jj^> a ^ are the modified widths, Qy , Q x are the momentum 
transfers, and F and A are the fragment and parent mass numbers. 

The first terms in (4.15a) and (4.15b) are the unmodified widths in the 
longitudinal and transverse directions given in terms of Fermi momentum 
as 


F(A-F) 

A 


2 

a 


o 



(4.16a) 


2 F(A-F) 2 

i "* l a 

x * o 


(4.16b) 


The extra factor of 2 in (4.16b) is due to summing over x and y 
directions . 

Moniz et al , 90 measured the Fermi momentum of a range of nuclei 
via electron scattering. We have developed in the last chapter the 
theory of longitudinal and transverse momentum transfer. The expres- 
sions for the longitudinal and transverse momentum transfer in nucleus 
nucleus collisions are given in equations (3.39) and (3.40). 
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Experimental data on the widths of longitudinal momentum distribu- 
tions of projectile fragments were measured by Greiner et al.^-9 The 
widths of the transverse momentum distributions , when resolved onto the 
detector plane were found to be similar to the longitudinal widths 
a(P u ) ~ a(P x ) - o( P y ) (4.17) 

within 10% suggesting isotropy of fragment production in the projectile 
frame . 

Brady et al.^6 measured the widths of transverse momentum distri- 
butions of projectile fragments in the reaction l||La (1.2 A GeV) + 

Projectile Fragment + X. These measurements were made at the Bevalac 
Heavy Ion Super-conducting Spectrometer (HISS) with the MUSIC detector. 
Fragment momentum widths o( P x ) were measured in the detector plane and 


fit to the expression 


a(P x ) 


F (A 1 F) (MeV/c) 
(A-l) 


(4.18) 


with variable a 0 . It was found that Goldhaber theory, based on inter- 
polated value of Fermi momentum Pp “ 250 MeV/c (- a 0 ) under- 

predicted the widths. Values of a Q « 169 MeV/c was necessary, im- 
plying an unreasonable value for the Fermi momentum of 377 MeV/c for 
l^La. Based on our theory of transverse momentum transfer, we show in 
Table 6 and Figure 27 a much improved agreement with the experiment. 
This shows that collisional momentum transfer is substantially respon- 
sible for increased transverse widths. 

These calculations were done using the following steps based on 
our expression for transverse momentum transfer. In Equation (3.39) we 
used a Woods-Saxon density for ^-^La and Harmonic-well density for 12 c 
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Table 6: Transverse momentum widths o'± (Qj_) in MeV/c of projectile 

fragments in the reaction 1: ^^La(1.2 AGeV) + ^ 2 C -*• Projectile Fragment 
+ X, calculated from equations (4.15b and (4.16b). The last column 
lists the modified widths. When squared, the third column yields the 
modifications due to momentum transfer. 



a(P x ) 

Goldhaber 

(MeV/c) 

112.0 

157 

.8 

192 

.6 

221 

.5 

246 

.8 

269 


289, 

.8 

308, 

6 

326 


342 


357. 

8 

372 


386. 

8 

398. 

8 

411 


423 


434 


445 


455 



a t (b) 
ours 
; MeV/c) 


15.6 

16.5 
18.0 

19.5 

21.0 

22.6 
24.36 
26.23 
28.2 

30.2 

32.6 

35.05 

37.6 

40.3 

43.15 

46.16 

49.3 

52.6 

56.04 


a(P x ) 
xpt fi 
(MeV/c) 
ref. 26 


*l(Q±) 








Table 6 (Continued) 


Projectile 
Fragment 
(mass no . ) 


*<*x> 

Goldhaber 

(MeV/c) 


a t (b) 

ours 

(MeV/c) 


o(P x ) 
expt fit 
(MeV/c) 
(ref. 26) 


119 

118 

117 

116 

115 

114 

113 

112 

111 

110 

109 

108 

107 

106 

105 

104 

103 

102 

101 

100 

99 

98 

97 


465 

475 

483.7 

492.5 
501 
509 
517 

524.3 

531.5 

538.5 
545 

551.6 
558 
564 

569.6 
575 
580 
586 

590.6 

595.4 
600 
604.3 

608.5 


59.60 

63.34 

67.03 
75.7 
79.78 
84.14 
88.46 
92.9 

97.3 
102 

106.4 
111 

115.6 

125.7 
119 

123.5 
127.9 

132.2 

136.7 
141 
145 

149.2 
153 


702 

716 

730 

743 

756 

768 

780 

791 

802 

812.5 
822 
832.4 
842 
851 

859.6 
868 
876 
884 
891 
898 
905 
912 
918 


*±(Q±) 

ours 

(MeV/c) 


660.3 

674.7 

687.3 

700.3 
713 

724.7 

736.5 

747.3 

757.6 

768.3 
778.0 
788 

797.5 

807.5 

814.3 

822.5 
830 
839 
846 
854 
861 . 

867.5 
874 


96 


612.5 


157 


924 


880.3 








Table 6 (Continued) 


Projectile 
Fragment 
(mass no.) 

a(P x ) 

Goldhaber 

(MeV/c) 

a t (b) 

ours 

(MeV/c) 

o(P x ) 
expt fit 
(MeV/c) 
(ref. 26) 

ours 

(MeV/c) 

95 

616 

160.5 

930 

886 

94 

620 

164.2 

936 

892 

93 

623.6 

167.5 

941 

897.6 

92 

627 

170.9 

946 

903 

91 

630 

173.8 

951 

908 

90 

633 

176.7 

955 

912 

89 

636 

179.5 

960 

917 

88 

638.7 

182.0 

964 

921 

87 

641 

184.5 

968 

925 

86 

643.6 

186.7 

971 

929 

85 

645.9 

188.7 

974.6 

932.7 

84 

648 

190.6 

978 

936 

83 

650 

192.2 

980.8 

939 

82 

651.8 

193.6 

983.5 

942 

81 

653.5 

195 

986 

944.5 

80 

655 

196 

988 

947 

79 

656.4 

197 

990 

949 

78 

657.6 

197.7 

992 

951 

77 

658.7 

198.3 

994 

952 

76 

659.7 

198.7 

995.5 

954 

75 

660.5 

198.9 

996.7 

955 

74 

661 

198.9 

997.7 

955 

73 

661.7 

198.8 

998.6 

956.6 

72 

662.2 

198.6 

999.2 

957.3 




Table 6 (Concluded) 


Projectile 
Fragment 
(mass no.) 

a(P x ) 

Goldhaber 

(MeV/c) 

a t (b) 

ours 

(MeV/c) 

a(P x ) 
expt fit 
(MeV/c) 
(ref. 26) 

<M_(Q±) 

ours 

(MeV/c) 

71 

662.5 

198.0 

999.6 

957.6 

70 

662.6 

197.6 

999.8 

957.6 

69 

662.6 

196.8 

999.8 

957.5 

68 

662.5 

196 

999.6 

957.2 

67 

662.2 

195 

999.2 

956.6 

66 

661.8 

193.8 

998.6 

955.8 

65 

661.2 

192.5 

997.7 

954.7 

64 

660.5 

191.15 

996.7 

953.4 

63 

659.7 

189.7 

995.5 

952 

62 

658.7 

188.0 

994 

950 







Transverse momentum widths, MeV/c 
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Figure 27: Transverse momentum widths of projectile fragments In the reaction 

,39 La (1.2 AGeV) + 1 *€ -+ Projectile Fragment ♦ X, where X Is 
unidentified. Harmonic-well density tvas used for 1J C, Woods-Saxon 
for Lanthanum. 
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where the parameters have been listed before in table 2. Root mean 
square momentum transfer Qj_ and squared rms momentum transfer Qj^ were 
calculated, and widths oj_(Qj_) were calculated using (4.15b). Since the 
experimental widths are in the detector plane and momentum transfer is 
calculated in the reaction plane, information on the azimuth of the 
reaction plane is necessary. Since such information is not available, 
we assumed that the momentum transfer in the reaction plane and detec- 
tor plane are identical. 


CHAPTER V 


CORRECTIONS 

Modifications to the above formalism will be addressed in this 
chapter. These are called (i) The Deceleration correction and (ii) The 
Coulomb correction. In theoretical calculations, the assumption of 
constant velocity is frequently made. Momentum transfer introduces an 
asymmetry into the problem; the assumption of constant velocity of beam 
nucleons is then strictly not valid. Our evaluation of momentum trans- 
fer in Chapter IV can only be correct if the corrections to our calcu- 
lations are small. As will become apparent, this correction, related 
to deceleration of the projectile due to momentum transfer is indeed 
small at high energies (but not necessarily so at lower bombarding 
energies) . Modification of the expressions derived previously will be 
made first. Numerical evaluation will be carried out for these cor- 
rections, labelled collectively as "The Deceleration corrections". 

Coulomb repulsion of the charges in heavy- ions will be treated 
next. Ever since the beginning of the science of heavy- ions, the 
interplay of nuclear and Coulomb effects has unveiled new insights. 

This is also true for the problem treated here. An analysis of elec- 
tromagnetic dissociation cross-sections at high energy reveals the 
importance of electromagnetic (EM) fields generated. The Weizsacker- 
Williams method of virtual quanta will be Used to show the importance 
of the EM fields. Momentum transfer calculations will be performed and 
compared with the nuclear contribution. Finally, modifications to the 


88 


C' ^ 



89 


point -Coulomb assumption will be made in favor of a more realistic 
charge distribution. 


V.a Deceleration Correction 

The assumption of constant velocity (per nucleon) is frequently 
made in intermediate and high energy nucleus -nucleus as well as hadron- 
nucleus collision calculations. The change in velocity is assumed 
small at these energies. Momentum transfer to the projectile, as 
evaluated in this work will enable us to evaluate the change in veloci- 
ty (per nucleon) thus testing the assumption of constant velocity. We 
shall develop, in addition a self consistent framework for evaluating 
momentum transfer. 

Consider a beam with kinetic energy per nucleon of T/A GeV. Then 
the velocity per constituent nucleon is 


with 


v * /3c - c 


(1 - l/ 7 2 ) 1/2 


T /A 

7*1 + ; M n - rest energy of nucleon 

M n 


(5.1) 

(5.2) 


Momentum transfer AP 


MeV 


c 


to the beam implies a change in velocity 


Av - A0c - per nucleon. (5.3) 

M n • A p 

We parametrize the changed velocity as follows 

v(z) - v + 5 (z) (5.4) 

where 6(z) is a path- dependent correction with 6(z) « v. 

Previous expressions for longitudinal and transverse momentum 
transfer read (with APm - Qn , AP^- Q x ) 
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AP|| - -(ApA^) Jd 3 £ p Pp(£p) J d3 ?T 


and 


^ J ImS (x'.fp.fT) 


dz' 


AP ± - -(A p A T ) J d 3 |^ /> p (^>) J d 3 |* T px(It) 

V, J Re t (x',fp,f T ) £il 


(5.5) 


(5.6) 


The bracketed expressions were z-integrable with a suitable choice 

the two-body amplitude t. Now because of (5.4), this is not so. 

For S( z) « v, the following expansion can be used 


dz _ dz dz 

v(z) v + 6(z) v 



(5.7) 


Equation (5.5) now can be written as 


AP|| - -CApAf) Jd 3 ?p Pp(lp) Jd 3 ?T PT(?T) 

m CO 

V I Im t (x'.fp.^T)— 

J -® v 

-V ?p f Im t (x' ,f p ,f T )fil . Hill + ... 1 
j -<o V V J 

(5.8) 

Note that we have obtained a series that takes into account corrections 
to our previous expressions, (5.5) and (5.6). The first term in (5.8) 
is' (5.5), as expected. Successive terms are corrections to (5.5) due 
to the asymmetry introduced as a result of momentum transfer. For S(z) 
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« v, the higher order corrections should be small at high 
and the series should converge rapidly. 

A reasonable ansatz for S( z) is 

— S (z) a V(b , z) 
dz 


The above equation (5.9) implies 

6( z) •» | v(b,z) dz 

«4 *(b,z) 

- * (b ’ z) . Av 

X(b) 


with Av as in (5.3) and the phase shift operator is 


X(b,z) - 


J. 


V(b,z)dz - - 


k J -<*> 


M 


V(b , z)dz 


rel 


Also note that 

S(z) 5(«) Av . A 

v ' v J - , as expected. 


v z 00 v 


Equation (5.8) can now be written as 


9 corr 

AP„ - AP |i - AP || +... 


Av 

where the correction, 0( — ) is 

v 


APjj° rr - -(A p A T ) J d 3 f p />p(|p) | d 3 |^ 


r X(b.z) dz' Av 


x(b) v v 


7 [ Im t (x' 

So J -m 


energies , 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 



92 

Numerical evaluation of the ratio appearing in (5.11), x(b,z)/*(b) 

has been carried out for various projectile target combinations for 

impact parameters ranging from 0-15 fm. The energy dependence of this 

ratio has also been checked. This ratio has been compared for + 

^Be (fig. 28) and 1^0 + 208p^ at incident energy of 2.1 A GeV. Two 

extreme values of impact parameter (b - 0,7 (fm) for ^Be and b - 0,11 

(fm) for 208pb) were chosen. The ratio is almost impact parameter 

independent for the Be target whereas differences of 15-20% are 

observed for the Pb target. Since the ratio , the above 

X(b) Av 

indicates that the z-dependent change in velocity 6(z) is more pro- 
nounced for a heavier target (208pb) than a lighter one (^Be) at the 
same impact parameter. This is in accord with our physical understand- 
ing since ^Be matter density is roughly Gaussian so that the ratio is 
independent of b. For a heavy target such as 208pb ( matter density is 
more appropriately of the Woods -Saxon type so that the ratio depends on 
the impact parameter. Thus the ansatz (5.9) is physically reasonable. 
The resulting asymmetry APjj orr was evaluated for ^0 (2.1a GeV) + ^Be 
and is plotted in figure 29 (the correction AP| orr has been multiplied 
been multiplied by 10 for display) . It can be seen that the correction 
at impact parameters 4.52 - 5.53 fm is merely 3% or less. At smaller 
impact parameters, however, the magnitude of the longitudinal momentum 
transfer as well as the correction to it are slightly greater, indi- 
cating the impact parameter dependence of these quantities. 

Another question that we can investigate is the energy dependence 
of the correction AP^ orr . The energy dependence of longitudinal 
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Figure 28: The ratio of phase-shift parameter x (b, *)/X (W *• • function of b and z 

(where z Is the distance along the beam direction and b Is the Impact 
parameter) In the reaction ’*6 (2.1 AGeV) *Be — ♦ Projectile Fragment + X. 



Correction (first order) to momentum transfi 



4.52 4.72 4.92 S.12 5.32 5.52 

Impact parameter, fm 

Figure 29: First order correction to calculated longitudinal momentum transfer as 

a function of impact parameter for 1 *0(2.1 A GeV)+*Be -♦ Projectile 
Fragment +X. The correction P^*®" has been multiplied by 10 for 
comparison purposes. 



Percentage change in velocity per nucleon 
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Impact parameter, fm 

Fig. 30 Ratio of change In velocity per nucleon to Incident velocity per 

nucleon (a* determined by the Incident beam energy) for '*0 (2.t A 
GeVH’Be — >Pro)ecti!e Fragment+X due to momentum transfer as a 
function of Impact parameter. 
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momentum transfer has previously been displayed in figure 25. It was 
observed that longitudinal momentum transfer saturates around -1 A GeV. 
At 0.2 A GeV, for example, Longitudinal Momentum Transfer (LMT) is -160 
MeV/c at 4.52 fm for + ^Be. At 3 fm, the correction AP| orr is -100 
MeV/c, which is -27% of the longitudinal momentum transfer of 370 
MeV/c. At 2.1 A GeV, there is only a -10% correction at 3 fm. 

Av 

Finally the ratio — is plotted as a percentage at 2.1 A GeV for 

v 

160 + ^Be. Note that Av has been defined in (5.3). The ratio is -1.6% 
at these impact parameters and higher at 3-4 fm. At 0.2 A GeV, this 
ratio is higher as expected, again indicating that momentum transfer 
introduces relatively greater change of velocity (per nucleon) at .2 
A GeV than at 2.1 A GeV . 

These effects should be experimentally observable at intermediate 
energy heavy ion collisions. Note that the higher order corrections in 
(5.7) have not been evaluated because their magnitude is expected to be 
small at high energies. At intermediate energy, these may not be small 
so that all the higher order terms may need to be taken into account. 

V.b Coulomb Corrections 

A complete treatment of the problem of momentum and energy trans- 
fer in relativistic heavy ion collisions must take into account Coulomb 
effects. In low-energy heavy ion collisions, Coulomb effects play a 
significant role. At relativistic energies, the importance of the 
electromagnetic fields generated can be understood by analyzing the 
Coulomb dissociation cross-sections in the fragmentation reactions^. 

It was found that these electromagnetic dissociation cross-sections are 


sometimes comparable to the nuclear contribution and may overwhelm the 
latter for high Z (where Z is the charge number) targets such as Lead 
or Uranium. This is specially true of single -nucleon knockout 
reactions^ 2 Q f the type Ap + Target ■+ (Ap-1) + X for projectiles such 
as 12 C(2 . 1 A GeV), 16 0(2.1 A GeV) and 12 C(1.05 A GeV) on targets 27 Al, 
6^Cu, lO^Ag an( j 208pb. Excitation of the giant resonances^ such as 
the Giant Dipole Resonance (GDR) and to some extent the Giant 
Quadrupole Resonance (GQR) and their decay contribute primarily to the 
one-nucleon knockout cross-sections, although the contribution of 
magnetic Ml resonance has also been pointed out. 93 The common proce- 
dure for the analysis of such Electromagnetic Dissociation cross- 
sections is due to Fermi, Weizsacker and Williams; it is known as the 
Weizsacker-Williams method of virtual quanta^. The EM fields gener- 
ated by a relativistic projectile are equivalent to two plane wave 
pulses of radiation P]^ and ?2 impinging on the target, P^ along the 
beam direction and ?2 transverse to it. The equivalent photon spectrum 
has been derived in many texts including Jackson; it is included in the 
Appendix for completeness. It can be seen that the photon number 
spectrum scales as Z 2 where Z^ is the charge number of the target; 
hence the overwhelming contribution for high Z targets can be under- 
stood. 

We shall treat Coulomb corrections at various levels beginning 
with the proton-proton interaction. This will be generalized, via the 
Weizsacker-Williams approach to heavy- ion collisions at relativistic 
energies . 
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Consider the two-body transition amplitude of this work for pp 


interaction. 
c pp (e,y) - 


This was given by 

-yfTa pp (e) [<* pp (e) + i] 
m 

- t(o) exp (-Ky 2 ) 



(5.17) 


with K - l/2Bp p (e) and t(o) is understood to be 

t(o) - yfl/m a pp (e) [a pp (e) + i] [2*B£e)f 3/2 (5.18) 


where a pp (e) , a pp (e) and B pp (e) are the pp cross-section, ratio of the 
real to the imaginary part of the forward scattering amplitude and the 
slope parameter respectively; y is the relative separation between the 
charges. With the parameters in Table 1 (with a pp( e ) - -.374), the 
transition amplitude at 2.1 A GeV is 

t pp (e , Y) = (35 - 94i) exp(-1.19 Y 2 ) MeV (5.19) 

Momentum transfer to the target proton is semiclass ically 


AP - 


Vt pp (Y 2 ) ; dz - vdt, v - 0c 

J -GO V 


r* r a 

J -«o 3b 


a ; + 


3z 


t(o) exp(-1.19 (Tj 2 + z 2 )) — 


<143 - 382 1) b e'< b2 ♦ z2>(l ' 19 > ; ™ 


(5.20) 


For I b I - 1 fm, z - 1 fm, 


AP = -13.2 


MeV 

c 


b + 35.4 


MeV 

c 


i b 


MeV 


- ( -13.2 b - 35.4 z) _ with ib - 

c 


-z 


(5.21) 
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MeV 

At a separation of -1.4 fm, an incoming proton imparts -13 of 

c 

transverse momentum due to strong interactions. The Coulomb repulsion 
of the protons contributes only 


AP - 


2e2 

bv 


MeV 


(5.22) 


of transverse momentum to the target. However, as the separation grows 
large, the nuclear and Coulomb contributions become comparable. The 
above calculation of the Coulomb effect was performed by evaluating 
the z- integrated force in an impulsive collision. Only the transverse 
electric field contributes, the longitudinal field's contribution 
vanishes due to symmetry. 

For heavy- ions, the above method for treating the Coulomb repul- 
sion can be generalized via the Weizsacker-Williams method. 

Consider the collision of heavy- ions with charge numbers 
respectively, with relative velocity v - /?c per nucleon. The EM 
fields generated at the projectile by the target (and vice-versa) can 
be found from equation (D.l) of Appendix D with q - Z^e. The momentum 
transfer in an impulsive collision is 

AP - Z p e [ * E tran (t) dt - 2 Z p Z T e 2 /bv b (5.23) 

• -00 


with the longitudinal contribution vanishing due to symmetry. 

The above equation is the basis of our calculations. For ^ C and 
16 0 projectiles (at 2.1 AGeV) on various targets ranging from Be 
through Pb, we have previously calculated the transverse and longitudi- 
nal momentum transfer in Chapter IV. Coulomb contribution is now shown 
in Figures 31-33 and compared with the nuclear contribution. As can be 



Transverse momentum tr 



Impact parameter, fm 

Figure 31: Transverse momentum transfer (MeV/c) to O Projectile by the ’Be 

target due to Coulomb Interaction compared with the nuclear 
contribution as a function of Impact parameter (fm) at A GeV. 
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Impact parameter, fm 

Figure 32: Transverse momentum transfer (MeV/c) to 1'rojccl I to by the r/ AI 

target due to Coulomb Interaction compared with the nuclear 
contribution as a function of Impact parameter (fm) at n (n< -2.1 A GeV. 
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10.66 10.86 11.06 11.26 11.46 11.66 


Impact parameter, fin 

Figure 33: Transverse momentum transfer (MeV/c) to Projectile by the Pb 

target due to Coulomb Interaction compared with the nuclear 
contribution as a function of impact parameter (fm) at E )n€ *2.1 A CeV. 
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readily seen for heavier targets, Coulomb effects are non-negligible 
and may modify the transverse momentum spectra of fragments substan- 
tially. Note that an impact parameter cutoff b m i n was included in the 
calculations. This was done to separate the strong interaction effects 
from the Coulomb. 

Inherent in the Weizsacker-Williams approach is the point-Coulomb 
assumption for the heavy- ions. For impact parameters b > b m £ n this is 
a reasonable assumption. However, for collisions with b < b m ^ n where 
bmin “ ^0.1 (Projectile) + Ro.i (Target) (i.e. the sum of 10% charge 
radii of the projectile and target), the point Coulomb assumption 
becomes questionable. This can be seen in the departure of the poten- 
tial for a unit charge in a uniform charge distribution from that of 
the point-Coulomb potential^ t.e. 


V(r) 



r < R 


r > R 


(5.24) 


with R the uniform charge radius of the distribution. This departure 
in the potential has been plotted in figure 34. The corresponding 
momentum transfer is reduced in magnitude as the overlap increases. 

For heavy ions described as two uniform charge distributions, extension 
of the above argument throws into doubt, for b < b m £ n , the Weizsacker- 
Williams approach based on the point-Coulomb field. 


Departure from the point Coulomb (Z.e.e.=l) 
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1 1 i » 1 * 1 , , , — i 

2 3 4 5 6 7 8 

Radial distance r, fin 


Departure of the Coulomb potential felt by an Incident charge within a 
uniform charge distribution ("•Pb, Z-82) of radius R-1.26(A**(1/3)) 
from the point Coulomb field. The constant C In V cou ,(r)-C/r has been 
set to unity. 


Figure 34: 




CHAPTER VI 


SUMMARY 

In this work, an optical model description of momentum (longitudi- 
nal and transverse) transfer and energy deposition in relativistic 
heavy- ion collisions has been presented within the multiple scattering 
theory framework. Longitudinal and transverse momentum transfer have 
been evaluated and compared with experiment for various projectile- 
target combinations. Momentum "downshifts" of projectile fragments in 
the collision of 12 C (2.1 A GeV) , l g 0 (2.1 AGeV) and 12 C (1.05 AGeV) 
with targets ranging from ^Be through 2 ®®Pb have been evaluated using 
the above theory and compared with target -averaged data from experi- 
ments. Transverse momentum widths of Lanthanum fragments in the reac- 
tion of l^La (1.2 AGeV) with 12 q have also been calculated. 
Quantitative as well as qualitative agreement has been found with 
experiment. Thus the theory of this work can account for many features 
of heavy ion momentum spectra in a comprehensive fashion. 

Energy transfer calculations can also be undertaken using the 
optical model description along the above lines. Such calculations are 
significant because energy transfer determines the various channels 
available for de-excitation (in the two stage Abrasion-Ablation model 
of relativistic heavy ion collisions). Sophisticated experiments in 
the future will hopefully address this question. Theoretical calcula- 
tions are already in progress in this respect and await comparisons 
with future experimental measurements. 
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We have provided, through "The Deceleration correction," a method 
for investigating and correcting the oft-used constant velocity assump- 
tion. Momentum transfer to the projectile could invalidate such an 
assumption. Based on our calculations, we find at relativistic speeds, 
that this change is minor. Specifically, we found that the first-order 
correction to our calculations were small at high energies (v “ c) . 
These corrections, however, were found to be substantial at lower 
bombarding energies, provided the impulse approximation remains valid 
at these energies. 

The very important question of Coulomb effects has also been 
addressed. This has been done within the Weizsacker-Williams method of 
virtual quanta. We find a substantial momentum transfer (transverse) 
due to repulsion of the charges. Use of more realistic charge distri- 
butions in place of the point-Coulomb assumption made in the 
Weizsacker-Williams approach resulted in reduction of the magnitude of 
momentum transfer as the collision impact parameters grew smaller. 

The future directions for research in this area remain open. Only 
the single scattering term of the multiple scattering series has been 
utilized so far. The importance of the double- scattering term and its 
physical meaning within the context of our work remains to be explored 
in detail. A comprehensive theory of energy -momentum transfer within 
the framework of energy-momentum conservation requires the relativistic 
theory of nucleus -nucleus interactions. Although a satisfactory theory 
of proton-nucleus multiple scattering formalism (Dirac phenomenology) 
already exists, similar approaches have had limited success for heavy- 
ions. An alternate approach may be the theory currently known as 
"Quantum Hadrodynamics"^ which is a relativistic field theory of 
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strongly interacting mesons and baryons . Semi-classical solutions for 
energy -momentum transfer can be obtained from the above theory for the 
collision of heavy- ions. Such approaches already show substantial 
promise. Future work should definitely proceed along these lines. 
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APPENDIX A 


DERIVATIONS (OF MOMENTUM AND ENERGY TRANSFER FORMULAE) 

In this Appendix, we shall go through the derivation of the final 
result for momentum transfer p^tran between the ath target constituent 
and jth projectile constituent (equation (3.31)), 


P«J . 
tran 


| d 3 ?p ,j(* p ) J dhj p 0 <? t) 


vf 


(A.l) 


and the total momentum transfer (equation (3.32)), 

-* a T A p . 

P -TV P aj 
tot L . . L . tran 
a-1 j-1 


(ApAx) J d3 Cp Pp(lp) Jd-^fx px(^t) 


j: 


V ^p I t ( x ^p»^T) 


dz 
v 

(A. 2) 


where the symbols have been explained in the text in Chapter 3 . 
Equation (A.l) was derived from the expression 


a t 

A p 




i 

A p 


p 

1 M 

X 

j-1 

tran 

- - < gp,o(£p) 

gT.o^T) 

- X . 

j-1 ? P-J 


r 00 

A X 

A p 






L 

X 

j-1 

X 

a-1 

c aj <*.«P.«T) 

dz 

V 

gp.o(£p)gT,o(£T)> 

(A. 3) 


where the projectile and target many-body wave functions • gp ^ anc * 
lgTo > are written in terms of single particle states in a Slater 
determinant as follows 
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^jl(^Pl) ^J 2 <£pi) •••• ^Ap(^Pl) 

^jl^P2) ^j2(^P2) ^JAp(fP2) 

Hitt? A ) £)2<£p A > ^jAp(^p A > 

(A. 4) 

and where the single particle states have been orthonormalized. The 
ground state single-particle densities are defined as 

p(l p) “ <gp,o(lp) 6 ^p * £pj> | gp,o(^p) > 




where £ means all permutations of £pj have to be taken. For a - 
Perm 

1,2... A t target and j - 1,2... A p projectile nucleons, momentum transfer 

due to collision of aj pair is P®^ and the total momentum is obtained 

c r<in 

by summing over all such pairs . 

Explicitly, we can write 



116 


^tran “ " ^ •*>«,.* I 

A p ! A x ! J o 


Perm 


'jl«Fl> ... 


*. A a a ) 

jAp s p p 


|d 3 | xl . . -d 3 | T , AT I 


Perm *kl ( ^T1 ) 


^kA T ( ^V 


"P _ A P A T « 

I V ip j I E J c a j (*>£ p ^ — 
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(A. 6) 
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— l l || « 3 f p j 
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*00 f 
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m-1 1 1 
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(A. 7) 


| d 3 lp Pj(lp) j d ^T P a ( ^ T) % J 


7 cp j.:-j ^-vv * 


(A. 1) 


— A n At — . 

P =. f i P a J 

tot u . tran 

J-1 a-1 


- - (A p A T ) J d 3 | p Pp(|p) | d 3 | T px(^T) J tCx'lp.lx) — 


(A. 2) 


as claimed. 

The derivation of the energy transfer in relativistic heavy- ion 
collisions is analogous to the above derivation for momentum transfer. 



We assume that the energy of the projectile can be evaluated by taking 
the expectation value of the sum of projectile energy operators between 
the full state vector (i.e. equation (3.14)) 

E - < (-) j H p | («) > ( A - 8 ) 

where Hp - Internal Hamiltonian of the projectile 


H - - JL f V 2 + Z Vij 

P 2m n j-1 i<j 


(A. 9) 


Explicitly, the above (A. 8) can be written as 

t 

j 


An Af 


E “ <gp,o(tp) gT.o(tl) I e ^ X expfi/v Z X f bjtj (x' , | p , ft) dz '^ 

I \ 1— 1 a— 1 J -co / 


— i V/ 1 Vi j 

2m n j-1 ? PJ i<j 


A? ^ 

exp (- i/v Z X f t K j (x 7 > £p > ?t) dz '\ 

\ j-1 a— 1 J -« / 


ik*x 

e 


gp,o(^p)8T,o(CT> 


(A. 10) 


Using the identity 
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yields 


exp[ iS] 


Ap __ 

— I 7 A + I nj 

2m n j-1 *Pj i<l 


exp[-iS] 


1 P _2 r 
— I + l V ij 

2m n j-1 ^PJ i<j 


i 

2m n 


S, l V 
J -1 ^Pj 


l2 

s, 

1 

CM 

a 

< 

w 

1 

2 


1 

1 

»-* 

^r> 

(_b 

1 


(A. 12) 


a x a, 


P r 00 


-III 

a-1 j-1 J - 


t aj (x',| p ,| T ) %- 


(A. 13) 


Using 

[ A , BC ] - [A,B]C + B[A,C] 

One obtains the operator expression 


(A. 14) 


1 

) 

1 


1 

I 

S -A 


s. 1 . 

j-i § Pj 



Similarly, 


"P 

+ l 

j-1 




Ap 


S, 

s. i n . 

j-i 5 pj 



— ^ 

_ 


(A. 15) 


can be evaluated using (A. 15). 
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Combining all the terms, one obtains 


E - <gp,o(lp)gT,o(l T ) H p gp,o(tp)gT,o(^T) > 


+ < gp,o(^p)gT,o^T) 

2m n 


+ _L_ <gp,o(lp)gT,o<?T> 


I . * s I gp,o<^p)gT,o<^T) > 

j -1 <PJ I 


2m, 


n 


im 


gp,o(tp)gT,o(^T)> 


+ j- <gp,o(^p)gT,o(iT) f ^ pj S * ^ pj | gp,o(*p>gT,o(*T» 


(A. 16) 


The first term in (A. 16) will be recognized as the initial energy 
of the projectile E 0) p ro j. This can be taken to the left hand side and 
energy transfer defined as 


E 


tran 


E 


E 

o 


(A. 17) 


The second and the fourth terms can be combined. The first term 
in the resulting expression can be converted into a surface integral by 
use of the divergence theorem. This will vanish if the single particle 
states vanish sufficiently fast at infinity. The remaining expression 
can be shown to contribute nothing to the excitation energy for even- 
even nuclei such as ^C, considered in this work. For odd-even 

nuclei, this contribution is 0(3_) , therefore this has been ignored. 

A p 

The final result for energy transfer is then 


E 


tran 


E 


E 

o 


“ ^gp,o^p)g^ 0 <CT) 

2m n 



gp,o^p)g^ 0 (Ct)> 

(A. 18) 


120 


The rest of the calculation is analogous to the momentum transfer 
calculations i.e. one obtains 


E 


tot 

tran 




“j 


tran 


(A. 19) 


where 

E tran ” Energy transfer in the collision of the a-j pair 

vv 

n 

(A. 20) 


I. 


2 


on J * xt 


and 


E tran“ < A p A T>f d ^P />p(lp) |d 3 ? T P T W 


u 


t (x',|p,| T ) — 

- oo v 


(A. 21) 


where pj (fp) and pa(£T) are the single particle density of the 
projectile and a*-* 1 target nucleon and Pp(£p), pt(?t) are c ^ e nuclear 
densities of the projectile and the target, as before- 


APPENDIX B 




V 

■P 


V 


and P 


'M 


f 


We want to discuss how our expression for momentum, equation 
(3.20) relates to Wilson's expression for potential, equation (2.41). 
We expect an analogy on physical grounds because momentum transfer in a 
collision is directly related to the gradient of the interaction in a 
one to one fashion. The coupled equations in (2.40) 


->2 ^2 

+ K ; 0 (x) - 

x m/i 


2M A A 

— — l V , , (x) , , (x) (B . 1) 

*-* mu m r n 1 ' ' T m' u r ' ' 


( A +At ) .>1 


m' fi r 


define a potential matrix V with matrix elements 

“ <g p,m^p* I V opt^ X ’ ^p’ ^ I g p,m' ^p )g T,/i' ^T ):> 

(B . 2) 


with 


V opt (x ’ V S T ) = X C aj<V x j> 


(B . 3) 


so that the matrix looks like (2.44) with 


V(x) - 


v oo,oo (x) 

v oi,oo (x) 

V 10,00 (x) 


V 00,01 (X) 

V 01,01 (X) 

V 10,01 (x) 


V 00,10 (x) 

V 01,10 (x) 

V 10,10 (x) 


(B . 4) 
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We introduced in (3.20) similar expressions for the momentum of 
the projectile and its matrix elements were defined in (3.20). 

Equation (3.20) defines a coupled problem because the approximate wave 
operator 


^(x, £ p . ? T ) “ exp 


if™ 

- 1 J.» V I v « T > dz ' 


(B.5) 


induces transition in the internal eigenstates g and s' of the 

6 P >m ^T.m 

projectile and target, respectively. Equation (3.20) defines a matrix 
P(x). 


P(x) 


P 00,00 (x) P 00,01 < ‘ X) 


P 01,00 (x) 

p io,oo (x) 


P 01,01 (X) 

P 10,01 (x) 


p oo,io (x) 

P 01,10 (x) 

P 10,10 (X) 


(B. 6) 


This expression for momentum 


^m/i.m'/i^ " <g p,m^p^ 


‘r . , . 

^1. V*’’ V «T> dz ' 


A P -*> 

•i I V 

j-1 ? P.j 


• -£ v i <*■• v 


x e 


| g p,m ,( V T ):> 


(B . 7) 


can be compared to Wilson's expression, (B.2) above. 
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One can recognize that there is a close connection between the two 

expressions. This connection should be one to one, since for each 

channel defined by the eigenstates m,m' of the projectile and of 

the target, an interaction V , # induces a momentum transfer 

P , , . 
m/i , m ' /i ' 

The Sudden Approximation 

The connection between our derivation and the more familiar 
formalism of time -dependent Schrodinger theory will be discussed 
briefly. 


Consider 



* s <t» 


H 


A A 

lMt)> - (H +H ) 
s o 1 


*s (t)> 


ft - 1 (B . 8) 


where ip (t)> is the time dependent Schrodinger state vector and where 

A A S 

H 0 and are the unperturbed Hamiltonian and the interaction, 
respectively. In the "Interaction Representation" the above 


reads^ » ^ ^ 


i — 

at s 


A A 

i H t -1 H t 

e ° (-H q +H 0 +H 1 ) e ° I ^(t)> 


- ^(t) ^(t)> 


with 


(B.9) 


iH t 

° I V c)> 


^ I (t)> - e 


(B. 10a) 
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a iH t a -iH t 

Hj(t) - e ° Hj* ° 


The time evolution operator is defined as 


(B.lOb) 


^ T (t)> - U(t,t ) (t )> 
i o 1 o 


(B.ll) 


Specifically, 


0j(co)> - U(®, -00) 


0 i (-oo)> 


(B.12) 


where U must satisfy 


i — U (t, t ) - H_(t) U (t, t ) 

at ° 1 


U (t t ) - 1 

O O 


(B . 13) 


Writing the above as an integral equation, one obtains 


U (t, t Q ) - 1 - £ J t Hj(t') U (t' , t Q ) dt' + 


(B . 14) 


On iteration, the above yields 


u (t, t o> 


. I (l\ 


n-1 


W- • • H I (t n ) 


(B . 15) 


The above is a formal solution to the time -development of the 
state vector from initial time t Q to final time t. It is a power 

A 

series in the strength of the interaction Hj(t) and it is time -ordered, 
i.e., earlier times occur to the left of later times. Explicitly, 


U(t, t Q ) - 1 




dt. 


O l— 


Hj^) H x (t 2 ) 


(B . 16) 
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Consider adiabatic switching of the interaction. Define^ 

A A 

H (t) - e' £ l c lH 1 (B. 17) 

which implies that instead of turning on the interaction at T - (t-t 0 ), 
turn it on and off slowly. As t -> ± ® 

A A A 

H x (t) - 0 and H - H q (B.18) 

and the solutions to the Schrodinger equation take the form 


A 

-iH t -IE t 

|lKt)> - e ° |* > - e ° |tf> (B. 19) 


with H 0 > - E ^ > for stationary states, 
o o 


Then in this limit 

A 

iH t 

*j(t) -e 0 ,*(t)> - |*> (B. 20) 


which is independent of time. Therefore, 


iV* 0 ! “ I *£ > 

|V* x (-»> | - | 4 >^> 


(B. 21) 


Next, the sudden approximation will be discussed to justify 
(3.19). We can write^ 
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A 

I ^1 (°°)> " U (°°i "°°) | (-*)> 




- £ 


e 



I 


n-1 


- oo 



Hj(t 1 ) Hj(t 2 ) . . . H^) 


(B . 22) 


It has been shown in reference [47] that a necessary condition for 
the sudden approximation to be valid is that the interactions at two 
different times commute 


H J (t 2 ) 


0 


(B. 23) 


so that time-ordering is unimportant in (B.22). Following reference 
[47] 


|^ I (®)> 


exp 


-i r H x (f) dt' 


|^ I (-®)> 


(B. 24) 


which is the desired result. 

Having derived (B.24) it is necessary to relate these equations 
(B.22) - (B. 24) to our previous discussion in equations in Chapter III. 
Essential to the derivation of equation (3.19) was the commutativity of 
interaction matrices at different points along a straight line 
trajectory 

0 



(B . 25) 
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In the Interaction Representation, this is entirely equivalent 


Hj ( t ' ) , H x (t") - 0 


(B. 26) 


provided one identifies dz' - vdt' and dz" - vdt" assuming constant 
velocity per nucleon. The latter can be justified on the basis of l 
high energy assumption. 

The physical basis for the sudden approximation in the high energy 
context is as follows; the collision time t co n « t nuc where t nuc is 
the period associated with orbital motion of the nucleons in the 
nucleus. This condition 


t ,, « t 
coll nuc 


(B.27) 


is met in relativistic heavy ion collisions. We shall argue this by 


calculating the ratio 


For incident beam energies above 1 


GeV/N the nucleon velocity « 0c - c. For the distance we shall 

take typical nuclear diameter « 6 fn. Then t ^ For 

c 

t nU c we will consider the period of nucleons in a Fermi gas i.e. 


t - (E/h)* 1 - l - t 


For a typical nucleus with Fermi energy 


40 MeV, the ratio is 


'coll _ 6 fm x 40 MeV 

: 2ir(ftc) 


fi c - 197.3 MeV.fm (B.28) 
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which is small enough to justify the Sudden Approximation. Although 
details of the nuclear diameter and Fermi energy may vary, at 
relativistic energies the ratio will be small. We have omitted the 
Lorentz factor in the above argument, but the conclusions remain 
unchanged because we have overestimated the Fermi energy (not all 
nucleons are at the surface) and have considered the diameter of the 
nucleus and not of the nucleons themselves. 



APPENDIX C 


ANALYTICAL RESULTS USING GAUSSIAN DENSITIES 
Analytical results are useful in getting physical insights as well 
as establishing the validity of numerical calculations. We present in 
this section analytical calculations of momentum and energy transfer 
using Gaussian densities. We know from Chapter III that momentum 
transfer between ath constituent of the projectile and jth constituent 
of the target is 


? aj - • J d3 «P | d3 £ i «T> V } 


fo J._ C «J <X ’ + ? P ' ?T> 


dz 1 


(C.1) 


where the notation was explained following equation (3.32). The total 
momentum transfer is obtained by summing over a and j as 


trans 


Z' Z P 

« j 




A p A T J d3 lp Pp(fp) J d3 ^T /’t(It) 


(C.2) 


v f i(x',rp.r T ) — 

J.oo K V 


Define a vector G - Gb b + G ± £ with (X, Y, Z) - (b, Z) (C.3) 


AAA A A 

_L 


with 


G 



f t(x' + £ p - 



(C.4) 
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so that equation (C.l) be rewritten as 

- - | d3f p Pj(fp) | d^ | ~ G (x' + ^ J 


P . 


and with 


(C . 5) 


t - - t(o) e 


■k(x' + £p - £T ) 2 


(C.6) 


where t(o) is assumed to be real. 


The components of <f can be written with £p >b (£ ± ) as the 

^ A A 

components of £ p along b (and perpendicular to b) and £x,b (C TJ _ ) 
the components of along b (and perpendicular to b) . 


- £<£> J* (2 k) (b + £ p>b - ? T .b) 

V 1 K 


- K 


(b + ^p,b ' ^T,b )2 *"^p,l 


and 


G, - 


t(o) 


£ (2*) « p .± - <T.X > 

-*<b + « p>b - « T , b > 2 -*« pii 


Define the thickness function 


-n 


Ft (fr> - I P T (fx,b + Z T ) dz T 


(C . 7) 


(C . 8) 


(C.9) 


where 


_p T (^ T ) - Aj. 


[tj] [■ rrl 


(C.10) 
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Fx(?t) “ a T- 


L 2w a^, 


exp 


+ 2 


Similarly 


Fp(fp) - A p . 


L 2 , 4J 


]• 


exp 


L 2 af J 


r 3£ 2 1 

P 


2 ai 


where a p , a^ are the root mean square radii of projectile and 
respectively . Therefore , 

P - - J d^fj Fj(^x) J d^£p Fp(£p) j^G (x' + £ p - Ct) 


or 


where 


p - - | d 2 f p Fp(fp) Q (b + e pb ) 


Q - J Fx(^t) j^ G ( x> + £p ‘ £t) J “ Q b ^ + Qj. 


with 


- J d ^T Fj(^x) G b + £p,b “ ^T,b) 


or 


Q b - 2 k A T Z ^LJL 
^ V i K 


2w a| 


i(t> + C Pl b; £p x ) 


(C.11) 

(C.12) 
target , 

J (C . 13) 

(C.14) 
(C . 15) 

(C . 16) 


where we have used (C.7), and 
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1<b+ < p , b ; f, ± ) - J | d ? Tb « T1 (b * c„ ib + « Ttb ) 


3 [4 + 4] “ < b + < pb ' <Tt> > 2 


2a T 

T 


V L* V * tJ 


(C 


On evaluation this yields (with a - 3/2a^) , 


I(b+ f p , b ; € pX ) - 


(a+K)' 


i b ^P.b| 


exp 


ax 


(a+K) 


{{ b + <p, b } 2 + 4x}] 


substituting (C.18) into (C.16) yields 


Qb - (2k) A T 

* v 


V K 


_J_ *_•_<)> + t 9th ) 

2 2 
2*a T (a+K) 


exp 




Similarly, using the same method as in (C. 16) - (C. 19) , we obtain 


Q ± - (2K) 


3 w a £ 


P.-L 


v ” Ac . 2 . , .2 

2?ra T (a+K) 


exp 


ax 


(a+K) 


{{ b ^p, b } 2 + «p,x}" 


• 17) 


(C.18) 


(C. 19) 


(C . 20) 


which yields the two components in (C.15). 
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Now we need to evaluate (C.14) 

- | d 2 ^, Fp(?p) Q(b + £ p ,b> - Pb b + P 


> Jp JT ^ 

> — y y p 

trans L ‘ u aj 

« j 


Pb b + P ± 1 


On evaluation this yields (a' - 3/2a^) , 


P b b - b A p A T (2K) 

V » K 


r^aa'b 


v ^ K 0 2 0 2 

2*a p 2*ax 


(aa' + k ( a+a ' ) ) 


or, 


exp 


aa r kb^ 


(aa' + K(a+a' ) ) 


P b b - b (A^) 2t(o) ^ (aa'b) x exp - xb 2_ J 


where x “ 


aa ' K 


aa ' + k ( a+a ' ) 


Thus, finally, 


T - 5 ; p - Pb b + p x 1 
trans 7 “j 

aj 


P b b and P x - 0 due to symmetry 


- b A p A T 2t(<>) aa'bx exp ^ - xb 2 J 


A 

1 


(C . 21) 


(C . 22) 


(C . 23) 


(C . 27) 


which is the analytical result for momentum transfer. 
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The above result holds for only real two-body amplitudes t(o) . For 
complex t(o), there is a longitudinal component of momentum transfer as 
we have shown in Chapter III. The above result can be generalized to 
the latter case. 


APPENDIX D 


THE EQUIVALENT PHOTON SPECTRUM 

In this Appendix, we shall derive via the Weizsacker- Williams 
method of virtual quanta the equivalent photon spectrum generated by 
the target at the projectile. 

For an incident particle of charge q, velocity V - 0 (c) passing 
a system S at an impact parameter b, the spectrum of equivalent 
radiation is obtained from the electromagnetic fields^ 

E 2 (t) - 5 * b 

(b 2 + 7 2 v 2 t 2)3/2 

E^t) - z (D.l) 

(b 2 + 7 2v2 t 2)3/2 

B 3 (t) - / 9 E 2 (t) ; ( 1 , 2 , 3 ) - (z,b, ^-b) 

where Ei(t), E2(t) are the electric fields along the beam direction and 
transverse to it and B3<t) is the accompanying magnetic field. For V 
~ c, the fields Ei(t), E2(t) and B3(t) are completely equivalent to 
plane wave pulses of radiation P]^ and P2 incident on S, P3 along beam 
direction, P2 transverse to it. The equivalency is not exact for P2, 
since there is no magnetic field accompanying E^(t). If the 

motion is nonrelativistic in the frame S, then the particles in this 
frame respond to electric forces only so that one can add an extra 
magnetic term to the fields without affecting the results. This field, 
B^(t) “ 0 E^(t) will be shown to be of minor importance in calcula- 
tions . 
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The equivalent frequency spectrum (energy per unit area per unit 


frequency interval) of the pulses P^, P 2 are given by 


dli(w,b) 

dw 2 it 

dl 2 ( w ,b) 


- f_ |E 2 (w) I 2 


(D.2) 


dw 


— |E!(w) 2 
2 n \ 


with 


E 2 (w) - — [ E 2 (t)e iWt dt 

J -“ 


(D.3) 


and similarly for E^(w) . Explicitly, these Fourier integrals are 


E 2 (w) 


qyb 


1 


iWt A*. 
e dt 


V7* (b 2 + 7 2 V 2 t 2 ) 3/2 


- - I 

r~ J 


® i wbx/yV , ... 

e ' ' dx yVt 

; x - ' 


>/T* - (1 


bV J -co 


+ x 2 ) 3 / 2 


(D.4) 


(D. 5) 


_q 

bV 


2 

7T 


7V l 7V J 


(D.6) 


and Ei(w) - - i 


_q_ 2 

7bV 7T 


— ^0 i — — \ 

7V \ 7V ) 


(D.7) 


where Kq(Ki) are the modified Bessel functions of zeroth (first) or 
der . 

The frequency spectra are 
dli(w,b) ! _2 ' - 2 • 2 


-If A k \ 


(D.8) 
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and 


dl 2 (w,b) 

dw 


Li! 
*2 c 





2 


4 



(D.9) 


These have been plotted in reference 94 (Figure 15.7). The 

intensity of pulse ?2 involves a factor of 1 and is of minor 

2 

7 

importance for V =* c. To obtain the energy incident per unit fre- 
quency interval, one sums the frequency spectra over all impact parame 
ters. This is 


dl (w) 
dw 



dl^(w) dl 2 (w) 

+ 

dw dw 


bdb 


(D. 10) 


where b m ^ n is a minimum impact parameter beyond which other interac- 
tions take over (strong interactions in heavy ion collisions, for 
example) . The result is 


dl(w) 


dw 



(O + i i 2 <o 


bdb; e 


wb 

7 V 



x Ko(x)Kx(x) - i ^ 2 x 2 (K^(x) - K 2 (x)) 


(D.ll) 


w ^min 
where x - 

7 V 


w « 



^min 


equation (0.11) reduces to 


For low frequencies 
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dl(w) _ 2 q 2 
dw it c 




(D.12) 


For high frequencies (w » 7 V/b m £ n ) , the result is an exponential fall- 
off 


dl 

dw 


i£l 2 /i . ivVc 2 ! 

w c \v t \ 2 ’ 


-2 


w ^min 

7 V 


(D- 13) 


The number spectrum of virtual quanta is 
dl (w) 


dw 


dw - N(ftw) d (ftv) ( ftw ) - E N(E)dE 


(D.14) 


so that 


N (B ) - 1 q2 


n (fie) 


q _ ( c ) 1_ 

(fie) \ Vf (fiw) 


x Kq(x)Ki(x) 


1 .2 2 , v 2 . , 

- P x (K. (x) 

2 


- Kp(x)) 


(D. 15) 


2 7 2 

— Z • a • 
?rE 




x Ko(x)K!(x) - - /9 2 x 2 (K 2 (x) - K 2 (x» 
2 i u 


(D . 16) 


wb 


For q - Z e, a - 


, x - 


min 


(ftc) 


7 V 


(D.17) 


From (D.9), the equivalent frequency spectrum for pulse ?2 

shown to contain a factor — . At relativistic speeds, this pulse ?2 

2 

y 

is, therefore, of negligible importance. This was the justification 
for adding a magnetic field Bi(t) - & Ei(t). The time integral of the 
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field Ei(t) was shown to yield zero by symmetry, so that this did not 
contribute to momentum transfer. The frequency spectrum confirms the 
above and justifies our retaining the transverse component in (5.21). 
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